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Abstract 

The Ramanujan sum Cn{k) and an{k), the kth coefficient of the nth cy- 
clotomic polynomial, are completely symmetric expressions in terms of 
primitive nth roots of unity. For k fixed we study the value distribution 
of Cn{k) (following A. Wintner) and an{k) (partly following H. Moller). In 
particular we disprove a 1970 conjecture of H. Moller on the average (over 
n) of an{k). We show that certain symmetric functions in primitive roots 
considered by the Dence brothers are related to the behaviour of Cp-i(fc) 
and ap-i{k) as p ranges over the primes and study their value distribution. 

This paper is an outgrowth of the M.Sc. thesis project of the sec- 
ond author, carried out under the supervision of the first author at the 
Korteweg-de Vries Institute (University of Amsterdam). Some of the nu- 
merical work described was done by Yves Gallot. 
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1 Introduction 



We recall from elementary number theory the notion of multiplicative order. If a 
and n are coprime natural numbers, then there is a smallest positive integer k such 
that = l(mod n), this integer k is the multiplicative order of a, ord„(a), and 
divides f{n), where is Euler's totient (discussed more extensively in the next 
section). We say that a is a primitive root of the multiplicative group (Z/nZ)* 
of units of Z/nZ if ord„(a) = (p{n). This means that (Z/nZ)* is cyclic and a is a 
generator. Let p be a prime (indeed, throughout we exclusively use the notation 
p and q for primes). It is well-known that (Z/pZ)* is cyclic and hence it has at 
least one primitive root. Indeed, if (7 is a primitive root mod p, then so are 
with j and p — 1 coprime and hence there are ip{p — 1) distinct primitive roots 
mod p. Let us denote them by gi, gt with t = ip{p — 1) and < gi < p — 1. In 
their paper j3] the Dence brothers consider symmetric functions of the primitive 
roots of primes. They consider for example 

S2{p) = Yl 3i9j (1) 

and show that this quantity, when considered mod p, assumes only the values 
in {—1,0, 1}. They write in particular (p. 79): 'As a matter of distribution, we 
observe that amongst the first 100 primes (beginning with p = 5) the residues 
—1,0, 1 of S2 occur in the ratios 12 : 59 : 29. The question of what these ratios 
should be in the limit of infinitely many primes is an interesting one'. The 
following table gives the ratios amongst the first lO-' primes with 2 < j < 6 (thus 
beginning with p = 2). 

Table 1: Value distribution of S2{p) 



71 (x) 


S2ip) = -1 


S2{P) = 


S2ip) = 1 




0.110000 


0.610000 


0.280000 


10^ 


0.099000 


0.625000 


0.276000 


10^ 


0.093000 


0.626100 


0.280900 


10^ 


0.094120 


0.627330 


0.278550 


10^ 


0.093939 


0.626216 


0.279845 



The sum in is the second elementary totally symmetric function of gi, ■ ■ ■ , gt- 
For > 1 let us consider more generally Sk{p) with Sfc(p) the kth order totally 
elementary symmetric function in gi, ■ ■ ■ , gt. 

Question 1. Fix k > 1. As p ranges over the primes, which values are assumed 
by Sk{p) {mod p) and with what frequency ? 

Likewise we consider, for A; > 1, the sum 

Sk{p) = E 9'- 
i<i<t 

For this quantity we consider the same question: 
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Question 2. Fix k > 1. As p ranges over the primes, which values are assumed 
by Sk{p) {mod p) and with what frequency ? 

One of the motivations of this M.Sc. thesis is to resolve these questions. In- 
terestingly, we will see that symmetric functions in primitive roots are closely 
related with some important objects in number theory; cyclotomic polynomials 
and Ramanujan sums. Both cyclotomic polynomials and Ramanujan sums arise 
in many contexts and hence we devote separate sections to them. 

We will see that the question of the Dence brothers naturally leads to the 
study of the value distribution of Ramanujan sums and of that of coefficients 
of cyclotomic polynomials. These two issues have not been well studied in the 
literature. Regarding the value distribution of Ramanujan sums we could only 
find an almost forgotten paper by the famous analyst Aurel Wintner [22] , which 
we reconsider in ^ Regarding the value distribution of coefficients of cycltomic 
polynomials, we reconsider a paper by Herbert MoUer ^21 in iHl and establish 
some new results. 

2 Preliminaries 

2.1 Multiplicative functions 

An arithmetic function is a function from the natural numbers to the complex 
numbers. A very important subclass of these functions are the so-called multi- 
plicative functions. These functions satisfy /(I) = 1 and 



The Mobius function and Euler totient function are important examples of mul- 
tiplicative functions. If / is a multiplicative function, then so is g{n) with 
di^) = J2d\nfi^)- Indeed, more generally, if / and g are two multiplicative 
functions, so is {f -k g){n) = J2d\nf(^)9(^/^)^ which is the so called Dirichlet 
convolution of / and g. If n = pl^ ■ ■ - pt" is the prime factorisation of n and / is 
multiplicative, then we have 



A closely related formula is the following one; if YldLi fi^) absolutely conver- 
gent and / is multiplicative, then 



f{mn) = f{n)f{m) for {m,n) = 1. 



s 



n(i+/(p)+---+/(p^o)- 



d\n 




The latter identity is known as Euler 's product identity. 
If / is multiplicative, then it is an easy observation that 



f{m)f{n) = f{{m,n))f{[m,n]) 



(3) 
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for all integers m and n, where (m, n) denotes the greatest common divisor of m 
and n and [m, n] the lowest common multiple. Likewise we define (mi, m2, . . . , m^) 
and [mi, 1712, . . . ,ms], were we put (mi) = 1 and [mi] = mi. Conversely, a func- 
tion that satisfies the functional equation Q is said to be semi-multiplicative. It 
can be shown [201 that if / is semi-multplicative, there exists a non-zero constant 
c, a positive integer a, and a multiphcative function fi such that 

^ |c/i(n/a) ifafn; 
1 otherwise. 

It is easy to see that if / is semi-multiplicative and A; is a constant, then also 
f{k/{k,n)) is semi-multiplicative in n. 



2.2 The Mobius function 

Let n be an integer having prime divisors pi, - ■ ■ ,ps- Then the Mobius fi function 
is defined as follows: 



a p'^\n for some prime p; 

(— 1)* otherwise. 



Note that fi{l) = 1 and that ^{nY = 1 iff n is squarefree. The Mobius function 
satisfies the following identity: 



E'-w-ii !,i:,L (4) 

d\n 



otherwise. 



An integer n is said to be kth powerfree if n is not divisible by a kth power of 
some integer > 1. Later we will need the following identity. 

Proposition 1 For any natural number k we have 



^/i(c/) 

|n 



1 if n is kth power-free; 
otherwise. 



First proof. Using the multiplicativity of the Mobius-f unction one sees that the 
left hand side is a multiplicative function of n. It thus suffices to evaluate it when 
n is a prime power. Now 



1 if a < A; - 1; 
otherwise, 



from which the result follows. □ 
Second proof. Write n = m^r, with the number r fcth power free. Then the sum 
under consideration equals Ylld\m^^^'^)■ '^^^ result then follows by (jD). □ 

An important formula involving the Mobius function is the celebrated Mobius 
inversion formula. It states that f{n) = J2d\n9i^) ~ ^d|n/^('^)/('^/'^)- 

The Mobius function often arises in combinatorial problems, we will see several 
examples in this M.Sc. thesis. In Section El some auxiliary functions involving 
the Mobius function are considered. 
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2.3 The Euler totient function 



The Mobius function often arises in combinatorial problems where so-called in- 
clusion and exclusion is being used in counting. We demonstrate this by deriving 
an identity for Euler's totient function ip{n), which is defined as 

An) = E 1- 

l<i<", U,n)=l 



Proposition 2 Let N{d) denote the number of integers 1 < m < n that are 
divisible by d. Then 

fin) = j:,(d)N(d) = "E ^ = "11 (l - ^) ■ 



q!|7i d\n p\n 



Proof. Suppose 1 < m < is an integer. In the expression '^d\nf^i^)-^i^)^ 
integer m is counted in those N{d) for which both d\m and d\n. It is counted 
with weight d\ml^i^)- -^^t' using (jH), we see that 



M = M = { 

',\n, d\m d\{n,m) 



1 if (n, m) = 1; 
otherwise 



Thus (f){n) = J2d\nf^(^)-^i^)- '^^^^ ^l^y then N{d) = n/d. On using that is 
a multiplicative function, the result then follows. □ 

From Proposition |21 we easily infer that is a multiplicative function. Alterna- 
tively this can be easily deduced on using the Chinese remainder theorem. 
For later use we note that by Proposition |21 we have 

V(n) ripiJl-i;) ,„ P 

where Si = Y[p\5 p\nP- Note that as n runs over all integers, 5i will run over all 
squarefree divisors of 5. In particular, (0) implies that 

(^(2n) r 2 if n is even; 
(^(n) 1 1 otherwise. 

2.4 Roots of unity 

An nth root of unity is a solution of 2;" = 1 in C. Note there are precisely n 
solutions; e^'^*/", ■ ■ ■ , e^™/", with, of course, e^™/" = 1. Instead of e^''*'^/" one 
often writes with ^„ = e^'^*/". Indeed, throughout this M.Sc. thesis we adopt 
the notation 

C,n = e . 

An nth root of unity is said to be a primitive root of unity if it is of the form 
with k and n coprime (which we will denote by (/c, n) = 1). Such a primitive root 
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of unity has the property that it does not satisfy an identity of the form = 1 
with m a divisor < n of n. 

Example. For n = 6, and C| = are the primitive roots of unity. Note that 
el is not a primitive 6th order root of unity as it satisfies the identity = 1. 

Note that there are (p{n) primitive nth roots of unity. The roots of unity form 
a cychc group of order n. The primitive roots of unity correspond with the 
generators of this group. 



3 Ramanujan sums 

The Ramanujan sum c„(m) is defined by 



l<fe<n l<fe<n 
{k.n) = l (fc,n) = l 



Although Ramanujan was not the first to work with Ramanujan sums, he was 
the first tor reahze their importance and use them consistently (especially in the 
theory of represenation of numbers as sum of squares, see e.g jHl Chapter IX]). 
The next proposition lists some of the important properties of Ramanujan sums. 
(These properties are all known or well-known except perhaps the formula in part 
6, for which we do not have a reference.) By t'p(n) we will denote the exponent 
of p in n, that is we have Upi^n) = r iff p^\n and p'^'^^ \ n. 



Proposition 3 

1) We have 



2) We have 



d\ (n,m) 



n ^ V'(^) 



Cn[m) = fi 

\ I 'T) m 1 / //-II 

^ (n,m) ' 

3) We have —{n,m) < Cn{m) < {n,m) and —^{n) < Cn{m) < f{n). 

4) We have Cn{m) G Z. 

5) We have c„^„2(m) = Cn,^(m)cn2(?T?.) if {ni,n2) = 1; i.e. c„(m) is multiplicative 
with respect to n. 

6) The function Cn{m) is semi-multiplicative in m. For fixed squarefree n, the 
function fi{n)cn{m) is multiplicative in m. For arbitrary natural numbers n and 
m we have 



7) The following orthogonality relations hold (when ri\r and r2\r): 



m=l 



(p{ri) ifri = r2. 
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Proof. 1) Consider g{n) := J2i<k<n^ """^ • Clearly this expression equals n if 
n\k and zero otherwise. Expressing g{n) in terms of Ramanujan sums we obtain 
9{n) = Y.d\n'^n/d{m), which by Mobius inversion gives c„(m) = Y.d\n9{d)Kn/d)- 
On using that g{d) = d if d\m and vanishes otherwise, the stated formula follows. 

2) This result, due to O. Holder, follows from property 1 on using (jSJ (see e.g. 
Exercise 1.1.14 of jl^). 

3) This is a consequence of property 2 and formula (0). 

4) First proof. Immediate by property 1. Second proof. The Ramanujan sum 
can be considered as an element of the field extension Q(Cn)- It is invariant under 
each of the automorphisms Cn "with k coprime to n. Hence it is in the fixed 
field of Q(Cn), which is Q, since Q(Cn) : Q is Galois. Since Cn is an algebraic 
integer, the Ramanujan sum must even be an integer. 

5) First proof. We have 

c„i(m)c„2(m) = 2^ e "i 2^ e "2 = 2^ e = c„,„2(m), 

l<fcl<ni I<fc2<n2 ^1,^2 

(fcl,nj) = l (k2,n2)-l 

where we use the observation that the set of congruences classes of the form 
kin2 + ^2^1, with 1 < kj < nj and {kj,nj) = 1 for j = 1 and j = 2, consists of 
ip{nin2) distinct congruence classes mod nin2, all of them coprime with nin2- 
Second proof. The assumptions on rii and n2 imply that {nin2, m) = [rii, m){n2, m) 
Now note that 

Cnin2M = ^ 5Z did2H{YJ^) = Cni{m)Cn2{m), 

d\{nin2,m) di\(ni,m) ^ 

d2l("2-'") 

where in the last step we use the multiplicativity of /i. 

6) The semi-multiplicativity is immediate from what has been said in Section ITT) 
and property 2. The second assertion is an easy consequence of property 2 and 
the observation that if / is multiplicative, then f{n)/f{n/{n, k)) is a multiplica- 
tive function in k. The proof of the third assertio uses in addition to the previous 
proof ingredients the multiplicativity of /i. 

7) For a proof see e.g. p. 17 of [22j. □ 

An arithmetic function / is said to be an even function of (n, r) if f{{n,r),r) = 
f{n,r) for all n (for a survey see 0]). By Property 2 Cn{m) = c„((n, m)) and 
hence Ramanujan sums are even. By Property 4 we can alternatively define 
Cn{m) = Ei<fc<n, (fc,n)=iCOs(27rm/c/n). 

An optimality property of Ramanujan sums was discovered by Bachman 
If r > 1 is any real number, he showed that for any sequence of real numbers bk 
we have 



I ^ ' - ^ 

m=l l<k<n \ / J 

(fe,n) = l 

It follows that if we consider the infimum of the left handside over all sequences 
bk with bk this is assumed in case bk = I for all k. 

For some further properties of Ramanujan sums, such as for example the 
Brauer-Rademacher identity, we refer to Chapter 2 of |l()j . 



n 

'cJm) 
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3.1 Density and order of growth 

is a set of natural numbers, by S[S) we denote the limit of x ^ X^n<a:, nes ^ 
(as X tends to infinity) if this exists. Let tt{x) denote the number of primes p < x. 
Recall that the prime number theorem asserts that asymptotically tt{x) ^ x/ logx 
(for a much stronger version see Lemma El below) . If 5* is a set of prime numbers, 
by 6{S) we denote the limit of 7i{x)~^ J2p<x pes ■'- (^^ ^ tends to infinity) if this 
exists. 

If / and g are two functions defined on a set 5", then we shall write f{x) = 
0{g{x)) (we say / is of order g) if the ratio \ f{x)/g{x)\ is bounded for all x E S. 
If the ratio f{x)/g{x) tends to zero for x tending to a specified value xq (which 
may be infinite), then we shall write f{x) = o{g{x)). The symbols O(-) and o(-) 
are commonly called Landau symbols. 



3.2 Ramanujan expansions of arithmetic functions 

Ramaujan sums play a key role in the theory of Ramanujan expansions. For 
completeness we discuss this topic, though this will be not needed in the rest of 
this thesis. 
Let 

M(/)= hm -J2fin). 

n<x 

denote the mean- value of / : N ^ C (if it exists). The Ramanujan sums have 
the following orthogonality property which is easily deduced from Property 7 of 
Proposition El 

ip{r) if r = s; 
otherwise. 

On certain spaces of arithmetic functions there is an inner-product {f ■ g) = 
M{f ■ g) and this product suggests Ramanujan expansions f ~ J2r>i^r{f)cr 
for arithmetic functions / in a suitable Hilbert space with coefficients ar{f) = 
(/ ■ Cr)/<^{r). Indeed, Ramanujan showed that for example 

(p{n) 6 jj,{r)cr{n) 



M{Cr 



n TT 

r> 



a.{n) = yd^ = ^y , MOc.(r^) (Re(,) > 1), 



and 



^in)=^^ -— : 



2r - 1 

where r(n) denotes the number of representations of n as a sum of two squares 
and C(s) as usual the Riemann zeta function (for Re(s) > 1 we have C,{s) = 
Yl'^=i and for all other s 7^ 1 in the complex plane ({s) can be uniquely 
defined by analytic continuation). Rearick [SOj proved that if F is an arithmetical 
function and an associated function / is defined by F{n) = Ylid\n^f^^) ^^'^ 
that Xln 1/(^^)1 < then we have, with a{n) = J2kf(^n), that the series 
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J2k 0'{k)ck{n) converges absolutely for each n and we have F{n) = J2k (^{^)ck{n). 
For a readable survey on Ramanujan expansions, see e.g. for more detailed 
information see [22] ■ 



3.3 Ramanujan sums and Dirichlet series 

Given an arithmetic function / an important quantity in analytic number theory 
is the associated generating series, Yl^=i f h^) I ^ which is usually named after 
Dirichlet. If we take /(n) = c„(m), we obtain 

oo , . 

n=l 

By property 1 of Proposition El we have c„(m) = ^^r=„ d\mf^i^)^ hence 



dr=n, d\m 



■s 



and so 



n=l dr=n, d\m d\m r=l ^ ^ 

(This proof is due to Estermann, for Ramanujan's (longer) proof see |H1 p. 140]). 
Taking s = 2 we obtain in particular 

71^ ( (-1)™ 2cos(2m7r/3) 2cos(m7r/4) 



ai(m) = ^ d = -m ( 1 + ^ + + 

which shows in a striking way the oscillation of cri(m) about its 'average' Tr^m/G. 
(It is not difficult to show that J2m<x "^i ("^) ~ T Sm<i ~ 12 ^^^'is to 

infinity: 

^.,(™) = ^<( = i^|£l(|£l + i) = i^H! + o(.^l), 

m<x rd<x r<x r<x \ r<x / 

for which it follows that J2m<x "^il"^) = 7r^x^/12 + O(xlogx).) 

Similarly one can consider -R„(s) = Ylm=i'^n{^)/^^ ■ Invoking property 1 of 
Proposition El again, one easily obtains that -R„(s) = ({s) J2d\n ^)^^~^ ■ 

In ^ the Dirichlet series fm{s) = ^"^=1 Cn{my / n"^ features (with j a non- 
negative integer). Define 



Using Properties 5 and 2 of Proposition IHl and Euler's product identity, one infers 
that, for Re(s) > 1, 



pS 
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It follows that 



= Fj^^\s)as)/a2s) and fj-\s) = Fl^^-'\s)/as). (7) 

4 Cyclotomic polynomials 

The nth cyclotomic polynomial, $„(a;), is defined by 

<l>n{X)= n iX-Ci). (8) 

l<j<n 

(j,")=l 

Note that the roots of $„(X) are precisely the primitive nth roots of unity. 
Clearly is a monic polynomial of degree (p{n). We write 

4>{n) 

<I>„(X) = ^a„(A;)x^ 

k=0 

It is the coefficients a„(A;) that later will have our special attention. Using in- 
duction with respect to n it is not difficult to show that an{k) G Z (see e.g. 
Exercise 1.5.24]). Note that 

^"-l = n n iX-0 = ll^n{X) = l[^,{X). 

d\n l<j<i d\n d\n 

(j,n) = d 

By applying Mobius inversion one infers from this that 

$„(x) = n(^^-ir^^^- (9) 

d\n 

Thus, for example, 

%{X) = = x'p-^ + + . . . + X + 1. 

Furthermore, it is known that ^2n{X) = $„(— X) if n > 1 is odd and $„(X) = 
$^(„)(X"/^(")), where j{n) = U, p is the squarefree kernel of n. The latter two 
properties are the reason that some authors in this area restrict themselves to 
n squarefree and odd. A further important property is that for n > 1 we have 
X'^(")$„(l/X) = <I>„(X). It is not difficult to infer these three properties from 
Q. In terms of the coefficients these properties imply (respectively): 

a„(/t) = ("^(")(^) if^l^; (10) 
1, otherwise, 

a2n{k) = {-l)''an{k) for n > 1, 2 f n; (11) 

(^n{k) = an{f{n) — k) for n > 1. 

It can be shown that $„(X) is irreducible over Q. This has as consequence that 
the degree of the cyclotomic field Q(Cn) over Q equals (^{n). Over finite fields 
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the cyclotomic fields factor in general. Indeed, over Fg, has 0((i)/ordg((i) 

distinct irreducible factors if (g, n) = 1. If p \ m, then if and only if the 

order of a(mod p) is m (Proposition ITHl below). This can be used to infer the 
infinitude of primes p = l(mod m) [16, Exercise 1.5.30]. 
Note that the coefficient, an{^{n) — 1), of X'^^'^'^^^ equals 

- ^ (n = -Cn(l) = 

l<j<ri 
(j,") = l 

where the last equality follows from Property 5 of Proposition El Since a„((y9(n) — 
1) = ctn(l) for ^ > 1, we conclude that an(l) = ~f^{^)- Another connection 
between cyclotomic polynomials and Ramanujan sums is given by the following 
result due to Nicol 

Proposition 4 We have 



oo 

r.Am \ 

m 



\ ^ — ^ m 

\ m=l 

and 

n 

J^c„(m)X™-i = (X"-l 



Using cyclotomic polynomials we can give an easy proof of Property 1 of Propo- 
sition ini We take the logarithmic of the right hand side of (jH)) and compare its 
Taylor series with that obtained on taking the logarithm of the right hand side 
of Q. The details are left to the reader. 



4.1 On the size of the coefficients of cyclotomic polyno- 
mials 

The size of the coefficients of cyclotomic polynomials is a much researched topic. 
This presumably stems from wonder over the fact that so often the coefficients 
are in {0, ±1}. Indeed, only for n > 105 some coefficients outside this range 
appear. For example, we have aio5(7) = —2. The amazement over the smallness 
of a„(m) was eloquently worded by D. Lehmer jH] who wrote: 'The smallness of 
a„(m) would appear to be one of the fundamental conspiracies of the primitive 
nth roots of unity. When one considers that a„(m) is a sum of C'^*^"'^) unit vectors 
(for example 73629072 in the case of n = 105 and m = 7) one realizes the extent 
of the cancellation that takes place'. 

Migotti showed in 1883 that apq{i) G {0, ±1}, with p and q odd primes. On 
the other hand Emma Lehmer showed that apqr{i) can be arbitrarily large (with 
p,q and r odd primes). Sister Marion Beiter put forward in 1971 the conjecture 
that in fact apqr{i) < (p + 1)/2 for all i and for any p < q < r and that this bound 
is best possible. This conjecture remains unsolved, for some recent progress see 

Let us put B{k) = max„>i |an(A;)|. For small values of k, MoUer [12] gave a 
method to compute B{k). The next table (taken from gives some values of 
B{k). 
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Table 2: Value of B{k) for 1 < A: < 30 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


B{k) 


1 


1 


1 


1 


1 


1 


2 


1 


1 


1 


2 


1 


2 


2 


2 


k 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


B{k) 


2 


3 


3 


3 


3 


3 


3 


4 


3 


3 


3 


3 


4 


4 


5 



The table suggests that B{k) > k for every A; > 1, this is however very far from 
the case: Bachman ITj extended work of several earlier authors and showed that 

logBW = C„,^fl+of™lV (12) 



(log A;) 1/4 \^ 

In the opposite direction, let us put A{n) = max^ |a„(m)|. Erdos has shown that 
there exists a c > and that there are infinitely many n such that 



log74(n) ^ exp 
On the other hand it is known that 



log A{n) < exp (log 2 + o(l)) 



clogn 
log log n ^ 

logn 



log log n 



where the constant log 2 is best possible. So the conspiration of the primitive 
roots of unity is not always that efficient ! 

Jiro Suzuki PH] proved the following beautiful result, which we present with a 
slightly shortened proof. Proposition lTTl using a slightly more involved argument, 
strengthens Suzuki's result. 

Proposition 5 We have {a„(fc) : n, G N} = Z. 

Proof. Given any integer s > 2 it is a consequence of the prime number theorem 
(cf. the proof of Proposition ^ that it is possible to find primes 2 < pi < p2 < 
■ ■ ■ < Ps such that pi + P2 > Ps- Let q be any prime > Ps- If s is even, let 
m = pi - ■ ■ PsQ, otherwise let m = pi ■ ■ ■ p^. We claim that 

$„(X) = J2 (^m{j)X' + (1 - s)XP^ (mod X^'^+i). (13) 

j=0 

The claim shows that am{p.s) = 1 — s. By (1111) we have a2m{Ps) = s — 1. On 
noting that 04(1) = the result follows. We next prove (fT^ . By and the 
observation that {d : d\m, d < ps + 1} = {l,pi, . . . ,Ps} and that for i 7^ j, 
Pi + Pj ^ Pi + P2 > Ps, we infer that mod XP'"*"! we have 

$™(X) = (l-XPi)---(l-XP=)/(l-X) (modXP=+^) 

= (1 + X + ■ ■ ■ + XP^-^){1 -X^' XP»-i) (mod XP»+^) 

Ps-l 

= ^"iij)X^ + (1 - s)XP» (mod XP^+^). 

j=0 

13 



This concludes the proof. 



□ 



The basic idea of this proof seems to have been first formulated by I. Schur in a 
letter to E. Landau. He used it to show that a„(A;) can be arbitrarily large. 

A very readable survey on properties of coefficients of cyclotomic polynomials 
is provided by 

In this thesis we are especially interested in how often certain values are 
assumed by a„(fc) for small fixed k (®, where n either runs over all integers n 
or over the numbers of the form p — 1 with p a prime ( ^9.3|) . 

5 Primitive roots 

We already stated that for every prime p there exists a primitive root g mod p 
and that all other primitive roots mod p are of the form with j coprime to 
p — 1. Since g, g'^, gP~^ are all distinct mod p, there are precisely ip{p — 1) 
distinct primitive roots mod p. 

Emil Art in conjectured in 1927 that for every integer g not equal to —1 or a 
square, there are infinitely many primes p such that g is a primitive root. More 
precisely, he conjectured that the set of primes p for which is a primitive root 
has a density which should equal a rational multiple of A, the Artin constant, 
which is defined as 

A = 17(1 - — -) = 0.3739558136 ■ ■ • 
p P{P - 1) 

Artin's conjecture was proved in 1967 under the assumption of the Generalized 
Riemann Hypothesis (GRH) by C. Hooley. In our consideration of the problem 
stated in the introduction, the constant A will also arise. Artin's constant (and 
indeed many similar ones) can be evaluated with many decimals of precision by 
expanding it in terms of values of the Riemann zeta- function evaluated at integers 
exceeding one, see e.g. ^3]. We return to the Artin primitive root conjecture in 
Section EZl 



6 Some sums involving the Mobius function 

Two important auxiliary quantities that will arise are, for given positive integers 
k, r and j, are the sums 

e^nd Yl (14) 

{n,r) — l {p—l,k) — r 

A number of the form p — t with t fixed is said to be a shifted prime. Note that we 
can restrict to the case where j = 1 and j = 2. The case j = 2 can be dealt with 
by fairly elementary analytic number theory. The case j = 1 is rather harder, in 
the case of shifted primes providing a non-trivial estimate would be considered a 
major achievement by the experts we consulted. 
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6.1 Counting squarefree integers 

In this section we are concerned with estimating the first sum appearing in ()14|1 . 



Lemma 1 Let r > 1 be an integer. We have 

6x 



7n<x J- ip|r V ' p> 

(m,r) — 1 



where the implied constant is absolute. 
Proof. We have, by inclusion and exclusion, 

(m,r) = l (d,r) = l 

where Ar{x) denotes the number of integers n < x that are coprime with r. Note 
that 

and hence Ar{x) = (f{r)x/r + 0{(p{r)). On using the latter estimate we obtain 
E M-)^ = x^Y: ^ + 0(v^^(r)). 



' r ^ d? 

m<x d<y/x 
(m,r) = l (d,r)^l 



X 



, E ^+0(/;».(r)). 

(d,r)=l 



0(v^^(r)), 



where we used that 

V{r) ^ Kd) <^{r) Yj,^ 1 , <^(r) 1 



^ ^ ^/ c(2)rn,,(i-^) c(2)n,,(i + j) 

and C(2) = vrVe. □ 
We also need to estimate Mr{x) = ^ ni<x jj,{m). Estimating this quantity turns 

(m,r) — 1 

out to be much harder. The quantity Mi{x) is the summatory function of the 
Mobius function and this function is called the Mertens function. The prime num- 
ber theorem is known to be equivalent with the assertion that lim^^oo Mi{x)/x = 
0. The celebrated Mertens' conjecture states that |Mi(x)| < ^/x for every x. If 
this conjecture would hold, then from the easily proved equality 

1 Ml (x) ^ 

-as. 



C(s) Ji X 



s+l 



it would follow that ({s) has no zeros with Re(s) > 1/2, i.e. the Riemann Hy- 
pothesis would follow. In fact, the Riemann Hypothesis is known to be equivalent 
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with the assertion that Mi(x) = 0{x^^'^^^), for every e > 0. In 1986 Odlyzko and 
te Riele disproved the Mertens conjecture. The behaviour of Mi{x) is thus closely 
related to the zero free region of the Riemann zeta function. The largest known 
zero free region of Ci^)^ due to Korobov, was used by Walfisz to show that 

Mi(x) = 0(xexp(-C(logx)^/^(loglogx)-^/^)). 

In this work we will be satisfied with merely determining the average of Mr{x): 

Lemma 2 We have 



m<.x 
( Tn , r ) — 1 



where the implied constant may depend on r. 

To prove the lemma, we will apply the Wiener-Ikehara Tauberian theorem in the 
following form. 

Theorem 1 Let f{s) = 'Yl^=i^nln'^ be a Dirichlet series. Suppose there exists 
a Dirichlet series F{s) = J2'^=ibn/n'^ with positive real coefficients such that 
(a) \an\ <bn for all n; 
(Jo) the series F{s) converges for Re{s) > 1; 

(c) the function F{s) can be extended to a meromorphic function in the region 
Re{s) > 1 having no poles except for a simple pole at s = 1. 

(d) the function /(s) can be extended to a meromorphic function in the region 
Re{s) > 1 having no poles except possibly for a simple pole at s = 1 with residue 
r. 

Then 

an = rx + o(x), X oo. 

n<x 

In particular, if f{s) is holomorphic at s = 1, then r = and J2n<x = o{x) as 
X — > oo. 

Proof of Lemma^ We apply the Wiener-Ikehara theorem with F{s) = ({s) and 



f( )= \^ ^ 1 



(n,r)=l 

Of course F{s) satisfies the required properties and has a simple pole at s = 1 
with residue one. Since the finite product in the formula for f{s) is regular for 
Re(s) > 0, the result follows on using the well-known fact that 1/C{s) can be 
extended to a meromorphic function in the region Re(s) > 1 (and hence r = 0). 
□ 

6.2 Counting squarefree shifted primes 

In this section we are concerned with estimating the second sum in 
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6.2.1 On a result of Mirsky 

Theorem |21 below is due to Mirsky ^21 (with the difference that in his version 
the 0-constant depends at most on k,r and H). In his paper Mirsky states 
two theorems, of which he only proves the first (the proof of the second being 
similar). Mirsky's second theorem is stated below. For completeness we give the 
proof here. Recall that the letters p and q are used to indicate primes. 

Theorem 2 Let r be any non-zero integer, k any integer greater than 1, and H 
any positive number. Then 

#{g < x : g - r is k-fiee} = 11 " pk~i^p-i) ) ^^^^^ + ^ (i^) ' 

where the 0-constant depends at most on k and H. 
Our proof rests on the Siegel-Walfisz theorem: 

Lemma 3 [19, Satz 4.8.3]. Let 7r(x;m, /) denote the number of primes q < x 
with q = /(mod m). Let C > be arbitrary. Then 

7r(x; m, I) = + 0(xe-'''^^), 

ip[m) 

uniformly for 1 <m < log*" a;, (Z,"^) = 1, where the constants depend at most on 
C. 

Note that if (/,m) > 1 there is at most one prime p = Z(mod m). If (Lm) = 1, 
then the above result implies that S{p = /(mod m)) = l/ip{m), thus asymptoti- 
cally the primes are equidistributed over the primitive congruence classes mod m. 

Proof of Theorem^ The dependence of 0-constants on r, k and H is indicated 
by using them as index, thus Oh means that the implied constant depends at 
most on H. 

Let y = \og^ X. By Proposition Q we have 

< a; : g — r is fc— free} = /i(a) 

a<y <i<x <^>y 

g^r (mod ) q = r (mod a^) 

= h + h, (15) 
say. We have, using Lemma 01 with C = kH, 

a<.y 

(a,r)^l g = r(mod a^) 



a<y 
(a,r) = l 
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E 4TYLi(^) + Oh{-^) + 0{y). (16) 



a<j; 
(a,r) = l 



Clearly 



a<y r-K J r-V / \a>v 

(a,r) = l (ii,'') = l 



By Euler's product identity we have 

/i(a) ^ TT / . 

(or^fc^ 111 



(a,r) = l ^1 

Using that /c > 2 and the classical estimate (p{m) ^ m/ log log m, we obtain 

On combining (fT^ with (fT7|) . (fTH|) and (fT^. we infer that 

■^f V p'' \p-l)J V (loga;)y / log^'^x 

Note that 

i^^isE E i£Ei^i = o(-^) = <). 

a>y m<x a>y 

m = r {mod a^) 

where in the second sum the summation is over all integers m < x satisfying 
m = r(mod a''). On adding the estimates for Ii and I/2I, the result follows from 

(uni). □ 



6.2.2 Connection with Artin's primitive root conjecture 

By Theorem |21 the density of primes p such that p — 1 is squarefree equals A, 
the Artin constant. The Artin constant also arose in Section in the context 
of primitive roots. This raises the question whether there is some connection 
between the two problems. We will now show that this is indeed the case; both 
problems are in fact special cases of a generalization of Artin's primitive root 
conjecture. We assume some familiarity with algebraic number theory. 

We recall that a prime p splits completely in Q(Cn) iff p = l(mod n). Note 
that p — 1 is squarefree iff p ^ l(mod g^) with q any prime. It follows that p — 1 is 
squarefree iff p does not spht completely in any of the fields Q(Cg2) with q a prime. 
It is a consequence of Chebotarev's density theorem that the density of primes 
q that split completely in a normal extension K : Q equals 1/[K : Q], where 
[K : Q] denotes the degree of the extension. On noting that the compositum 
of Q(Q), ...,Q(Cg|), where qi,...,qs are distinct primes, equals Q(C{gig2-g.)2)5 we 
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expect by inclusion and exclusion that the density of primes p for which p — 1 is 
squarefree equals 

h [Q(0) : Q] ^ ^(^') 

Indeed, the primes p for which a given integer g with (7 7^ —1 or a square is 
a primitive root mod p can be described in a similar way. Here we want that 
for each prime q with g|p — 1 we have that g 1 ^ l(mod]9), which can be 
reformulated as p does not split completely in Q{Cq^9^^'') for S'liy prime q. By 
inclusion and exclusion we expect then that that the density of primes p for which 
g is a primitive root equals 



n=l 



^^"^ ^- (20) 



In fact Hooley proved in 1967 that the number given in ()2Up is indeed the 
correct density, under the assumption of GRH. For squarefree n the degree 
^n, g^^^) '■ Q] is generically equal to nLp{n), which is the exact degree of 

and hence the arisal of A in both problems does not surprise us. 
Both problems considered here are in fact special cases of the following gener- 
alization of Artin's primitive root problem. This generalization was first studied 
by Goldstein and later by M. Ram Murty ^Tj. Let K be an algebraic number 
field. Let JF be a family of number fields, normal and of finite degree over K. 
Determine the number of prime ideals V of K such that Nk/q{V) < x and which 
do not split completely in any element K of J-'. On taking K = Q and F to be 
the set of fields of the form Q{Cqi 9^^'')^ where q runs over the primes, we obtain 
Artin's primitive root problem. On taking K = Q and F to be the set of fields 
of the form Q(Cg2), where q runs over the primes, we obtain the special case of 
Mirsky's result considered above. 

From an heuristical point of view the arisal of Artin's constant in the square- 
free problem is not surprising. The density of primes q such that q ^ l(mod p"^) 
is 1 — l/{p{p — 1)). Imposing the condition q ^ l(mod p"^) for each prime p and 
assuming all the conditions to represent independent events, one arrives at A as 
the expected density of the primes q with q — 1 squarefree. The independence of 
the various local conditions is a consequence of the fact that for any two distinct 
primes r and s we have Q(Cr2) H Q(Cs2) = Q- 

For the Artin problem we impose the local condition that q does not satisfy 
both q = l(mod p) and g^'^'^^P = l(mod p). The density of primes q not satis- 
fying the latter two conditions is 1 — l/[Q(Cp, fi'^'^^) : Q]- If the conditions for the 
various primes p would be independent, we would expect (as Artin original did) 
the density to equal 

1 



n 



(21) 



which if (7 7^ —1, 0, 1 equals a rational number times the Artin constant. However, 
it is not always true (with r and s as before) that Q{Cr,9^^^) H Q{Cs,9^^'^) = Q 
(for example when r = 5, s = 2 and g = 5 in which case the intersection 
equals Q("\/5)) and hence we expect the density to come out different from (j2T| 
sometimes (as can be shown true assuming GRH). 
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6.2.3 A variation of Mirsky's result 

We will need a variation of Mirsky's result for k = 2 and r = 1. Let S = 
(gi, . . . , qt) be a (possibly empty) finite sequence of primes satisfying qi < q2 < 
... < qt. We write I'si.'n) = (ug^i^ri), . . . ,h'q^{n)). We say that a number n is 
S-squarefree if f tt, for every prime q with q ^ S. We put 



(— 1^^^ if n is S-squarefree; 
otherwise. 



If S is the empty set, then Hsin) = fi{n). 

Proposition 6 Let t he any non-zero integer, Let S = {qi, . . . ,qt} and ci, . . . ,et 
be natural numbers. Write Q = 115=1 ^Ij^ ■ have 



p<a; 



n n 



where the O -constant depends at most on S,ei, . . . ,et and H. 

Proof. A variation of Theorem |21 and left to the reader. □ 

Remark. Again the constant appearing in this result is not surprising. For each 
prime q we impose a condition at our prime p occurring in the sum: if g G 5 we 
prescribe z/g(p — 1) and if g ^ S" we require that p ^ l(mod g^). If for fixed g we 
compute the densities of primes satisfying this condition and if we multiply all 
these 'local densities' together, we arrive at the density 

t 

6{f^s{p-l)j^0)l[6{u,^{p-l)=e,) = 



n(--^)n(-^)n 



e,=0 



which equals the density given in Proposition IHl 
6.2.4 A conjecture 

It is generally believed that primes p are such that p—1 does not have a preference 
with regards to having an even or odd number of prime factors. In other words, 
the following conjecture (thought to be deep by the experts) is generally believed: 
we have J2p<xf^iP ~ ^) — x tends to infinity. We propose a slightly 

more general conjecture. 
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Conjecture 1 Let t be any non-zero integer, Let S = {qi, . . . ,qt} and ci, . . . ,et 
be natural numbers. We have 

1) = o(7r(x)), 

'^S(p-l) = (<!l,....et) 

where the o-constant depends at most on S,ei, . . . ,et and H. 

All the numerical data we came across in this respect seemed to be not inconsis- 
tent with this conjecture. 

7 Value distribution of Ramanujan sums 

We discuss the value distribution of c„(m) as n varies over the integers and m 
is fixed. We do this following Aurel Wintner's paper |2S1- Wintner argues that 
the prime number theorem is equivalent with the Mobius function having an 
asymptotic distribution function and notes that Cn(l) = /i(n). A natural question 
arising then is whether c„(m) for m fixed has an asymptotic distribution function 
and if so, what it looks like. We next recall some facts regarding moments and 
distribution functions. 

7.1 Distribution functions 

A sequence {aj}°lo ^^^^ have an asymptotic distribution function, p, if there 
exists a monotone function p{a) such that 

< 1 

lim — — = p(a), 

holds at every continuity point a of p and, moreover, we have p(— oo) = and 
p(oo) = 1. It is known that if limj_,oo sup \aj\ < oo, it has an asymptotic distri- 
bution p iff Mn{aD exists for every j; in which case 

/oo 
a^dp{a). 
-oo 

(We define a° = 1 and hence [p]o = 1.) Given a distribution function p(a) we 
can consider, L{u; p), the Fourier-Stieltjes transform of p: 

/oo 
e*""dp(a). 
-oo 

If the moments exist for every j, then 

is valid for every u. In the case where there are only finitely many distinct 
elements in ths sequence (which will happen in our application), p(a) will be 
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a step function making only finitely many steps. Then L{u] p) will be a finite 
expression of the form Yl^j=i c(/5j)e*^'' and the values assumed are then Pi, - ■ ■ , P^, 
each with density c{(3j). 

As an illustration let us compute the asymptotic distribution for {fi^n)}'^^^. 
By now the reader should have no difficulties in proving directly that this is 
merely the step function which has the saltus ('jump') S/vr^, 1 — Q / tc'^ , 3 / n"^ at 
a = —1,0,1, respectively. Let us derive this by the method of moments. The 
odd moments of n are zero and the 2jth moment with j > 1 equals G/vr^. The 
zeroth moment equals one. We thus derive that 

r/ \ 1 6 6 3 ^ 6 3 ,„ 

L{u,p) = 1 - ^ + — cosM = — e + 1 - — + — e*", 

from which we arrive at the same conclusion as before. 



7.2 Computing the moments and asymptotic distribution 
function 

Using the representation (|7j) for the Dirichlet series /P^(s) and f!^^^~^\s) and the 
Wiener Ikehara theorem. Theorem^ the mean of the jth moment of c„(m) can be 
determined (note that by property 5 of Proposition El we can take F{s) = m^({s) 
in that result when we apply it to fn\s)), We find that 

Mnicnimfn = 4^ = ^ n^,^XU(l)- (22) 

Note that M„(c„(m)^-')(^(2) is a multiplicative function in m. Similarly, we have 
Mn{cn{m)'^^~^) = 0. (We use the subscript n in M to indicate that the mean is 
with respect to n.) Since |c„(m)p-' < rn^^ and hence the same inequality holds 
for the mean, it follows that the radius of convergence of the power series 

^ M„(c„(m)^^') 2 . 
Gm{u) = 2^ 7— i-u y (23) 

j=o ^ •'>■ 

is infinite for every m (= 1,2, ■ ■ ■). In case m = p"^ , Gm{u) is easily evaluated 
using and (jHl). One obtains 



tt'^ I ~ p J ' n'^ \ p J \ p^ p^'^^ 



Note that if k is fixed, then 



6 6 

Gpr{u) ^ ^ COSM, (p ^ oo) = G'i(m), 



and 



cos(^(p^>_) ^ (A: ^ oo). 



.=0 ^ 
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For any given m it is not difficult to compute Gm{u) using formula (j22|) . The 
existence of all the moments M„(c„(m)-'), j = 0, 1, 2, ■ ■ ■ implies the existence of 
an asymptotic distribution function p^- Note that L{u; pm) = Gmiu) and that 
Gm{u) = X]j=i c(/3j)e*^''. The values assumed are Pi,---,Pr, each with density 
c{(5j). Thus, taking m = for example we immediately read off from the formula 



for Gpr{u) that 



if^ = 0; 



5{cn{f) =v) = l ^2pfe(Vi) if ^ = ±</'(/), where /i = 1, . . . ,r; 



2 , 1, if t; = ±p''. 



8 Value distribution of cyclotomic coefficients 

8.1 Evaluating cyclotomic coefficients 

Write 

We consider the value distribution of the cyclotomic coefficients a„(A;) for fixed 
k, as n runs over all the positive integers. By setting fi{n/d) = whenever n/d 
is not an integer we have, for n > 1, 

00 

$„(X) = - X'^)'^^"/'^) = JJ(1 - (24) 

d\n d=l 

Note that, for IXI < 1, we have 



00 00 / ^ 00 

_ = n 1 - M^)^'^ + ^M^)(M^) - 1) E 

d=l d=l \ j=l 



(25) 



where we used the observation that, for |X| < 1, 

{l-Xr^'^/'^=(^l-p{^)XU^-p{^)ip{^)-l^^^ (26) 

From (j25|) it is not difficult to derive a formula for a„(/c) for a fixed k; this is just 
the coefficient of in the right hand side of (j24j) (this approach seems to be 
due to D.H. Lehmer P). We thus obtain, for n > 1, 

a„(l) = -/i(n); 

a„(2) = fi{n){^{n) - l)/2 - fi{n/2)- 
a„(3) = yu(n)V2 - yu(n)/2 + /i(n/2)p(n) - /i(n/3). 

More generally, we have 

an{k) = J2 <ku fc.; ei, . . . , es)/i(^)'^ ■ ■ ■ /^(^)'% (27) 

where the sum is over all partitions ki + ■ ■ ■ + kg of all the integers < k with 
ki > k2 > ■ ■ ■ kg and over all ei, . . . , with 1 < ej < 2 for 1 < j < s. Using this 
result we will deduce that the n dependence of a„(/c) is not that strong. 
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Proposition 7 Put = \cm{l, 2, k)Y[p<kP- can uniquely decompose 
n as n = UkCk with {ck,Nk) = 1 and and Ck natural numbers. There exist 
functions Ai and Bi with as domain the divisors of Nk such that 

^ / Ai{nk)iJ,{ckf + Bi{nk)iJ,{ck) if nk\Nk; 
\ otherwise. 

Proof. The assertion regarding the uniqueness of the decomposition n = rikCk is 
trivial. For a given n only those partitions ki, k2, ■ ■ ■ , kg contribute to (j27j) for 
which n/ki is an integer for 1 < i < s. Note that ki\nk. Thus, we can write 

^^f)" ■ ■ -^^f)" = ■ ■■Kfri^ickr^-'-''- 

rC]^ rCg nJ]^ rCg 

If Uk \ Nk, then none of the integers Uk/ki, ...,nk/ks is squarefree and a„(A;) = 0, 
so assume that nk\Nk. On using that fi{r)'^ with w > 1 either equals /x(r) or 
fi{rY, the result follows from (jTTj) . □ 

The above proposition shows that {an{k)\n G N} is a finite set, thus if we fix 
k, there are only finitely many possibilities for the values of the coefficient of 
in a cyclotomic polynomial. We will now show that —1,0 and 1 are always 
amongst these values. In the formula for a„(A;) there is always the term —^{n/k). 
Let us take n = c/cj^^^^p, where c only has prime divisor > k. Then all the 
terms of the form ix{n/r) with 1 < r < are zero (since either r \ n oi n/r is 
not squarefree) and we obtain that a„(A;) = — /i(c)(— 1)'"^'^^ where 7r(x) as usual 
denotes the number of primes p < x not exceeding x. In particular, it follows 
that an{k) always assumes the values —1,0 and 1. Combining this insight with 
the previous proposition we arrive at the following result. 

Proposition 8 Let k > 1 be fixed. Let Nk = \cm{l, 2, k)Y[p<kP '^'^^ 
q > k be any prime exceeding k. We have 

{-1, 0, 1} C {0, a,(fc), ad,{k) \ d\Nk} = {a,(A;) | n > 1}. 

Using formula ()10|) it is seen that in the latter proposition A^^ can be replaced by 
kY[p<kP- We have |a„(A;)| < max„>i |an(^)| = B{k). See Table 2 for the values 
of B{k) for 1 < A; < 30. 

We next show that the inclusion in Proposition |S1 is strict for k > 13. Our 
proof rests on the following rather elementary result on prime numbers. 

Proposition 9 For k > 13 there are consecutive odd primes pi < P2 < Ps such 
that ps < k < pi + p2. 

Proof. Breusch proved in 1934 that for a; > 48 there is at least one prime in 
[x, 9x/8] (this strengthens Bertrand's Postulate asserting that there is always a 
prime between x and 2x, provided x > 2). Let a = 1.32. A little computation 
shows that the above result implies that for a; > 9 there is at least one prime 
in [x, ax]. One checks that the assertion is true for k G [13,21). Assume that 
k > 21 {> 9a^). Let ps be the largest prime not exceeding k and let pi and p2 be 
primes such that pi,P2 and p^ are consecutive primes. Then p^ > k/a, p2 > k/a^ 
and pi > k/a^. On noting that Pi + P2 > k{l/a + 1/a^) > k, the proof is then 
completed. □ 
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Proposition 10 For k > 13 we have {—2, —1, 0, 1} G {a„(A;) : n G N} (and thus 
B{k) > 2). 

Proof. Let Pi,P2 and be odd primes satisfying the condition of Proposition El 
Using we infer that 

where we computed modulo X'^"*'^. It follows that ap^p^p^{k) = —2. The proof is 
completed on invoking Proposition |H| □ 



8.1.1 Numerical evaluation of an{k) for small k 

For our purposes it is relevant to be able to numerically evaluate a„(/c) for small 
k and large n. A computer package like Maple evaluates a„(A;) be evaluating the 
whole polynomial For large n this is far too costly. Instead it is more 

efficient to use ()24p and expand for every d for which fi{n/d) 7^ 0, (1 - 
as a Taylor series up to 0{x^^^) and multiply all these series together. The most 
efficient method to date is due to Grytczuk and Tropak f^. First they apply 
formula (HUJ. Thus it is enough to compute a„(A;) with n squarefree. If < k, 
then an{k) = 0, so we may assume that 4>{n) > k. Let d = {n,Y[p<kP)- 
Tr = fi{n)fi{{r, d))ip{{r,d)). Compute bo, . . . ,bk recursively by 60 = 1 and 



1 ' 

bj = — : ^ bmTj-m for i > 0. 



? 

m=0 

Then 6^ = a„(/c). The proof uses the formula 



o-nik) = — — an{m)cn{k — m) for k > 1, (28) 

which follows by Viete's and Newton's formulae from (jH)) and it uses Property 2 
of Proposition El An alternative proof of ()28|) is obtained on using the following 
observation together with Property 1 of Proposition |21 

Proposition 11 Suppose that, as formal power series, 

00 00 

(1=1 d=0 

then dr{d) = X;^=i r{d - j) Y.k\j 

Proof. Taking the logarithmic derivative of ndlill ~ X'^)~°-'' we obtain 
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whence the result follows. 



□ 



For every integer f it is a consequence of Proposition^lthat there exists a minimal 
integer k, k^in, such that there exists a natural number n with an{kmin) = v. 
Grytczuk and Tropak [71 Table 2.1] used their method to determine fcmin for 
the integers in the interval [—9,..., 10]. Gallot has extended this range from 
[—60, . . . , 70]. From Table 11 one can determine fcmin for the range [—15, ... , 15]. 

8.2 Moller's paper reconsidered 

In this section we reconsider Moller's paper |13j. In Moller's approach a„(/c) is 
connected with partitions of k. A partition of a positive integer m is an expression 
of the form mi + ^2 + ■ ■ ■ + = m with all the rrij > 0. Ordering is disregarded. 
Thus, 1 + 1 + 2 + 3 and 3 + 2 + 1 + 1 are considered to be the same partitions 
of 7. A partition can be identified with a sequence {nj}'^Q of non- negative 
integers satisfying = W.l.o.g. we can denote a partition. A, of k as 

A = (fci 1 where > > . ■ ■ > > 1 (thus the number kj occurs 

times in the partition). The set of all partitions of m will be denoted by 
V{m). The number of different partitions of m is denoted by p{m). Hardy and 
Littlewood in 1918, and Uspensky independently in 1920, proved that 

p{m) ~ e^v/^/(4mV3) as m ^ oo. 



Proposition 12 ^13j Satz 2]. We have 



Proof. The taylor series of (1 — X)" equals, for \X\ < 1, Yl^oi~^)'' i'j)-^'' ^ where 
(") = a(a — 1) ■ ■ ■ (a — (j — l))/j!. Using this we infer that 

(1 - = X^(-iy (^^^^^ |X| < 1, (29) 

j=o \ 3 J 

The proof now follows from ()29p and ()24|) . □ 



Our proof above shows that the formula given in Proposition El is a triviality, 
whereas Moller's ingenious and rather involved proof of it obscures this. 
MoUer uses Proposition to show that 

T\J ( tlW r ^n<x (^n{k) 

Mn{an{k)) = hm = 

exists and gives a formula for it. To do so he first determines the average of 
He does this by expanding it out as a sum of terms of the form 
(n/j) with 1 < ctj < 2, which is possible by (|Hn|) . The average of each 



26 



term Yl fJ'"'^ [n / j) is easily determined. He then invokes Proposition IT2l and uses 
combinatorics to simplify his expressions. His final result is still quite complicated 
and we did not feel inclined to check its equivalence with our more simple formula 
in Theorem El below. We will now follow MoUer's trail, in spirit, but not in detail 
and see where it leads us. 

Proposition 13 Let s > 1, ki, ■ ■ ■ , kg be distinct integers and > > • • • > 
^ 1- Iff^ki >2 we let t be the largest integer < s for which rik^ > 2, otherwise 
we let t = 0. Let L = [ki, . . . ,ks] and G = {ki, . . . ,ks). We have 

lim 1 = 

' \nkj \ ) vr^ Gn,|^(p+1) ' 

where 

{1 if '^'ki = 1, s is even and fi{L/G) ^ 0; 

fi{^y^y2 ifuk, > 2 and niL/ki) = ■■■ = fi{L/kt) and ^{L/G) ^ 0; 
otherwise. 

Corollary 1 We have 

n<x ^ ' ^ ' 





Proof of Proposition El Put 

.w.E(-ir''--(1f)---(1f). 

n<x ^ ' ^ ' 

Comparison of and yields 

1 ifj = 0; 

-li{n/d) if i = 1; 

^ ^ ( I^{n/d){i2{n/d) -l)/2 ifj>2. 

Note that for j > 2 the binomial coefficient is only non-zero if fi{n/d) = — 1. 
Using (jHUj) it follows that a necessery condition for the argument of S{x) to be 
non-zero is that L\n. Now write n = mL. Note that fi{mL/ki) ■ ■ ■ n{mL/ks) = 
fi{mY fi{L / ki) ■ ■ ■ fi{L / kg) if {m,L/kj) = 1 for 1 < j < s and equals zero oth- 
erwise. It is not difficult to show that [^, . . . , -^j = ^ and using this, that 
fi{L/ki), ■ ■ ■ , fi{L/ks) are squarefree iff L/G is squarefree. It follows that if 
ji{L/G) = 0, then S{x) = and we are done, so next assume that ii{L/G) ^ 0. 
We infer that 



iy[^^'^^\ = {-Kn/d) ifj = l; (30) 



m<x/L, (m,L/G)=l 
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Let us first consider the (easy) case where > 1. Then we obtain S{x) = 
(-l)X^) ■ ■ ■ /^(^ ) ^m<x/L, (m,L/G)=i /^M'- If « is odd, then by Lemma H it 
follows that lim^-^oo S{x)/x = and we are done, so next assume that s is even. 
Then we apply Lemma ^ and obtain that 



The assumption fi{L/G) 7^ implies that + 1/p) = GY[p\l/g(p ~^ -'-)• 

Next we consider the case where n^^ > 1. The corresponding binomial coefficient 
is only non-zero if fi{mL/ki) = — 1. Similarly, we must have fi{mL/kj) = —1 for 
1 < j < t. It follows that if it is not true that fi{L/ki) = . . . = fi{L/kt), then 
S{x) = and hence limx^oo S{x)/x = as asserted, so assume that = 
. . . = fi{L/kt). We have 

s{x) = Yl (-M^)rv(^)---M^) 

m<x/L, {m,L/G) = l ^ ^ 

fi{mL / } — — 1 

^ ^^-^ * m,<x/L, (m,L/G) = l 

p{m) = -p{L/fc]^) 

On invoking Lemma and Lemma |21 the proof is then completed. □ 

Proof of Corollary ^ Follows from a case by case analysis from Proposition 
on using (j30|) . □ 

Theorem 3 We have 



where 



<A) = n(-i)% (^l^^)+f[{-ir^ ( 



L(A) = [/ci, ...,ks] and G{X) = {ki, ...,ks). 

Proof. The result follows from Proposition together with Corollary □ 



In case is a prime, G{X) = 1 for every partition and the above formula further 
simplifies. The above formula suggests a connection with the group or represen- 
tation theory of the symmetric group 5*^. The conjugacy classes in Sk are in 1-1 
correspondence with the partitions of k. If A = {k^''^ . . . k^'"' , then the order of 
every element in the corresponding conjugacy class equals L{X). In particular, 
L{X) < g{n), where g{n) denotes the maximum of all orders of elements in 5*^. It 
was shown by E. Landau in 1903 that log g{k) ~ y/ k log k as k tends to infinity 
(for a nice account of this see whereas by Stirling's theorem log k\ ^ k log k. 
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From Proposition 1121 MoUer infers that |a„(fc)| < p{k) — p{k — 2). To see this 
note that the partitions having 1 occurring at least twice do not contribute if 
/i(n) G {0,1}. If /i(n) = —1, then either /i(2n) = 1 or /i(n/2) = 1. This in 
combination with (fTTj) allows us then to argue as before and leads us to the same 
bound. By a much more involved argument, using some analytic number theory, 
Moller concludes that B{k) > for k > ko{m). This result was sharpened by 
several later authors culminating in Bachman's estimate (jl2j) . 

Moller used Proposition to evaluate B{k) for 1 < k < 30. He did this 
by hand. The outcome is in Table 2. We redid this computation by computer 
and arrived the same result. Note that M„(a„(A;)) = 6efc/7r^, with a rational 
number. For 1 < < 20 we give the value of Cfe in Table 3 (our table agrees 
with the one given in J^, except for the incorrect values eio = 319/1440 and 
ei6 = 733/2016 appearing there). 



Table 3: Scaled average, = C(2)M„(a„(A;)), of a„(/c) 



j 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 







1 


1 


1 


1 


7 


1 


7 


19 


31 




2 


6 


3 


8 


24 


18 


24 


144 


160 


j 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 




1 


55 


13 


61 


2287 


733 


667 


79 


55 


221 




16 


192 


288 


288 


20160 


4032 


8064 


336 


1344 


960 



Regarding Moller proposed: 

Conjecture 2 [13J. Let k>l. Write Mn{an{k)) 

1) We have < Cfc < 1/2. 

2) We have {-lf{ek - Ck+i) > 0. 



Moller stated that with help of an IBM 7090 he wanted to check his conjecture 
for further values of k. Had he carried this out, he would have discovered that 
(—1)^^(634 — 635) = —18059/4626720 < 0. Other counterexamples occur at A; = 
35,45 and 94. Indeed, we would not be surprised if part 2 of the Conjecture is 
violated for infinitely many k. 

On the other hand, part 1 of the Conjecture is true for k < 75. The numbers 
Cfc seem to be decreasing to zero and their size seems to be related to the number 
of prime factors of k, the more prime factors the larger seems to be. In Table 
11 (kindly computed by Yves Gallot) we give values for for k up to 61. 

8.3 Average and value distribution 

We give, using Proposition [7| a simpler formula for M„(a„(A;)) involving a„(A;) 
for a finite set of n. 



Theorem 4 Let k > 1 be fixed. Put = k Y[p<kP^ '^'^'^ let q > k be any prime. 
Then 

3 -sp ad{k) + adq{k) 

'1 + i) ^ 



M„(a„(fc)) 



p<k \ 



d 
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Furthermore, when v ^ 0, 



llp<fcV d\Mf. d\Mf, 



Proof. Let A^^ = lcm(l, 2, ■ ■ ■ , fc) np<A:P' '^i = Ilp<kP ^^'^ ^2 = riNk. We have 



n<x d\Nk 

(n.r2) = d 



J2 iMd)l^irnY + Bi{d)fx{m)) 



d|Affe m<x/d 

{m,r2/d) = l 



(m,r-|^ ) — 1 



d\Nk m<x/d 
(m,r| ) — 1 

where we used Proposition [7| and LemmaEl On invoking Lemma^we then obtain 
that 

E/7 \ 6 Ai(d) . , 

anik) = E ^ + ^^(-)- 

On noting that Ai{d) = {ad{k) + adq{k))/2 (and Bi{d) = {ad{k) — adq{k))/2, but 
this is not needed), the first formula follows, but with A^^ instead of M^. On 
invoking formula (fTUI) it is easily seen that 

^ k ~ ^ k 

d\Nk d\AIk 

and hence the first formula follows. The proof of the second formula is similar 
and left to the reader. □ 

Using identity pip we arrive at the following corollary to this theorem: 
Corollary 2 Let k > 3 be fixed and odd and Mk and q as in Theorem^ Then 

ad{k) + adq{k) 



2<p<k\^ ' P' d\Mk/2 ^ 



Furthermore, when v ^ 0, 



Sia„ik)=v)= ' ( E ^+ E i+ E ^+ S 25'- 

ilp<A;V py d|Mj./2 <i|A/j,/2 d\Mi./2 d|Mj./2 

As a corollary of Theorem 0] and Corollary El we have: 

Proposition 14 We have ek'^kY[p<k{P +1) G In case k is odd we have 
ekkUp<kiP+^)^^- 
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Yves Gallot observed that actually for k < 100, we have Ckk np<fc(P + 1) G ^• 
In case k is prime, the divisor sum in the previous corollary can be further reduced: 

Proposition 15 Let k > 3 be a fixed prime. Put = Y[2<p<kP^ ^'^^ let q > k 
be any prime. Then 

ad{k) + adq{k) 



MJaJk)) = — — ^ V 



2<p<fcV-^ ' p' d\Nk ^ 

Proof. We consider the formula given in the previous corollary. The divisors of 
Mfc/2 are either of the form d with d\kRk or of the form dk"^ with d\Rk. For 
the latter d we find, using (fTUj) that a^k-^ik) = adk{l) = —fi{dk) and hence 
^d\Mj2i(^d{k) + adq{k))/d = J2d\kRS^d{k) + adq{k))/d. Now suppose that d\Rk. 
Using Proposition El we infer that 

adk{k) + adkq{k) = ad{k) + adq{k) + ^i{d) + ^i{dq) = ad{k) + adq{k). 

Using this observation it follows that 

E ad(fc) + ddqik) _ ,\sr^ <^d{k) + adq{k) 
d ~^^k^^ d 

d\kRk d\Rk 

whence the result follows. □ 

Clearly 6{an{k) = 0) = 1 — X^^^o '^('^"(^) ~ '^)' where the latter sum has only 
finitely many non-zero values and it is a finite computation to determine those 
V for which an{k) = v for some n. For 1 < < 16 the non-zero values of 
({2)6{an{k) = v) are given in Table 4 (except for v = 0). A more extensive 
version of Table 4 is provided by Table 11. 



Table 4: Value of C(2)5(a„(A;) = v) 





v = -2 


v = -1 


V = 1 


v = 2 


k = 1 





1/2 


1/2 





k = 2 





1/12 


7/12 





k = 3 





5/24 


3/8 





k = A 





1/6 


1/2 





k = h 





13/80 


23/80 





k = Q 





25/144 


67/144 





k = 7 


1/576 


577/2688 


731/2688 


1/1152 


k = 8 





1/8 


5/12 





k = 9 





65/384 


347/1152 





k=10 





161/960 


347/960 





k=ll 


1/2304 


8299/50688 


11489/50688 


1/4608 


k=12 





349/2304 


1009/2304 





k = 13 


43/48384 


219269/1257984 


277171/1257984 


43/96768 


k = U 


13/21504 


2395/21504 


2319/7168 


1/2304 


k = 15 


13/32256 


1345/7168 


97247/322560 


13/64512 


k=16 


5/21504 


12149/64512 


1127/3072 


5/2688 
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Let us now look at Theorem IHl and Theorem 0] from the viewpoint of computa- 
tional complexity. In TheoremlHlthe sum has p{k) terms and the famous estimate 
of Hardy and Littlewood yields that log p{k) ~ ir^/TkJS as k tends to infinity. 
In Theorem^we sum over t{k) terms where logt(A;) ~ 7T{k) log2 ~ /clog2/ logfc. 
So Theorem El yields the computational superior method. Theorem |3] is, however, 
much more easily implemented. Using Proposition 1121 in combination with Theo- 
rem 0] a result comparable with Theorem |21 is obtained. Indeed, if one starts with 
Proposition and invokes Proposition ^[ one obtains a sum over partitions of 
k, where now 2 for example does not occur in the partition. This yields a result 
superior in complexity to that provided by Theorem 01 

As already pointed out by Moller one could use his method to study the 
value distribution of a„(A;) in case e.g. B{k) = 2 by considering the integer 
an{k){an{k) — l)/2 to determine 5(a„(A;) = —1) for example. This then yields a 
sum with p{k)'^ terms and this results in an algorithm that has worse complexity 
than that provided by Theorem |3J Aside from this, this seems to be, from the 
practical point of view, an unwieldy method. 

8.3.1 Some observations related to Table 4 

In this section we make some observations regarding Table 4 and Table 11, which 
is a much more extensive version of Table 4. 

Let us put B{k) = {a„(/c) : n G N}. Recall that \B{k)\ = B{k). Recall also 
that —1,0, 1 G B{k). Using this and Corollary El we infer that if k is odd, then 
B{k) is symmetric: we have v G B{k) iff —v G B{k). For k is even numerical 
results suggest that often B{k) is not symmetric. For A; < 75 it is true that if 
V G B{k) and v is negative, then —v G B{k). An other observation that can be 
made is that for < 75 it is true that B{k) consists of consecutive integers, i.e. 
if fo < vi are in B{k), then so are all integers between vq and vi. 

Let us define Bj{k) = {a„(A;) : n = j(mod 2)}, for 1 < j < 2. It is easy to see 
that always G Bj{k). Note that for k is even we have Boik) C Bi{k) = B{k). 
For k is odd we have v G Bi{k) iff —v G i3o(fc). We could also express this as 
Bo{k) = {—V : V G Bi{k)}. Thus in this case, if we know Bo{k), we also know 
Bi{k). Inspection of Table 11 shows that for odd integers k with B{k) > 2 often 
6{an{k) = B{k)) and 6{an{k) = —B{k)) differ by a factor two. Regarding this 
situation we have the following result: 

Proposition 16 Let k be odd and v ^ 0. We have 

26{an{k) =v) = 6{an{k) = -v) iS v ^ Bi{k) (that is iff -v^ Bo{k)). 

Proof. An easy consequence of Corollary |21 '<^=' This uses in addition that 
{anik) : 2 t n} = {adik) : d\Mk/2}. □ 

Example. Inspection of Table 4 shows that the condition of the proposition is 
satisfied for k = 7 and v = 2. It thus follows that there is no even integer n 
for which a„(7) = —2 (whereas aio5(7) = —2). Further examples can be derived 
from Table 5. 
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Table 5: Set theoretic difference B{k)\Bo{k) 



k=7 


{-2} 


A; = 11 


{-2} 


k=13 


{-2} 


= 15 


{-2} 


k=17 


{-3} 


A; = 19 


{-3} 


k = 21 


{-3} 


A; = 23 


{-4,-3} 


k = 25 


{-3} 


A; = 31 


{-4} 


= 35 


{5} 


A; = 37 


{5} 


A; = 39 


{5,6} 


A; = 43 


{-7} 


= 45 


{-7} 


A; = 47 


{-9,-8} 


= 51 


{8} 


A; = 53 


{9,10,11,12,13} 



For Proposition^Jto be of some mathematical value we would hope that infinitely 
often Bo{k) is strictly contained in B{k). The next result shows that a2n{k) 
assumes all values as n and k range over the integers. 

Proposition 17 We have {a2n{k) : n, k E N} = Z. 

Proof. Given any integer s > 2 it is a consequence of the prime number theorem 
that it is possible to find primes 2 < pi < pi + 2 < p2 < ■ ■ ■ < Ps < Pi + P2 — '2'- 
Let k = pi + P2 and let q be any prime > k. If s is even, let m = 2pi ■ ■ -psQ, 
otherwise let m = 2pi ■ ■ - ps- We claim that 

= ^a„(j)X^- + {s- l)X^-2 - {s - 1)X'-' (mod X'). (31) 

j=0 

The claim together with 04(1) = yields the result. We next prove 1)311) . Using 
the observation that {d : (i|m, d < k} = {1, 2,pi, ■ ■ ■ ,ps, 2pi} we infer from 
that mod X'^ we have 

= ^ : (mod X") 

(1 + X)(1-Xp2)---(1-XpO 
= (1 - X + X^ - X^ ■ ■ ^(l + XPi + ■ ■ ■ + X^^) (mod X'^) 

= Y,aUj)X' + is- 1)X'-' - {s - 1)X'-' (mod X'^). 

This concludes the proof. □ 
The above argument shows for example that 

'^'2-19-23-29-31-37-39(-^) = ■ ■ ■ — 5X^° + 5X^^ + • • • 
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9 Value distribution of Sm{p) and Sm{p) 

9.1 Ramanujan sums, cyclotomic polynomials and primi- 
tive roots 

Proposition El connects the mod p reductions of totally symmetric functions in 
primitive roots to totally symmetric functions in primitive roots of unity. Our 
proof rests on the following simple result. 

Proposition 18 Letp he a prime not dividing m. Then p\(^m{(i) iff 'the order of 
a(mod p) is m. 

Proof. If the order of a(mod p) equals m, then a"' ^ l(mod p) for every d < m 
and a™ = l(mod p). By © we then infer that Since there are ip{m) 

elements of order m in (Z/pZ)*, the equation = O(modp) has at least 

ip{m) distinct solutions mod p, however since the degree of equals (/^(m), 

these are all the solutions. □ 

Remark. We think this proof is more enlighting than that given in [16,, p. 209]. 
A related, less well-known result is that a{aP~^ — l)\E'p_^(a)/\E'p^i(a) is congruent 
to — l(mod p) if a is a primitive root mod p and divisible by p otherwise |18j . 

Proposition 19 Let p he a prime. Put t = ip{p — 1). Let gi,---,gt he the 
modulu p distinct primitive roots. Let 1 < ji < j2 < • ■ ■ < Jt < P — 1 he the 
natural numhers coprime to p — 1. Let f{yi,---,yt) he any totally symmetric 
function in the variahles yi, - ■ ■ ,yt. Then 

/(^/i,---,^7t) = /(C^i,---,C-i) (modp). 

Proof. By Newton's result on totally symmetric polynomial functions it is suf- 
ficient to prove the result for all elementary totally symmetric polynomial func- 
tions. By Proposition ITHl we infer that 

^p-i{x) = {x- gi)---{x- gt) (mod p). (32) 

On the other hand we have 

%^,{x) = {x-CiU)---{x-CiU)- (33) 
From (jH^ and (jHSJ, the result easily follows. □ 

Let gi,---, g^p(p-i) denote the distinct primitive roots mod p with 1 < gj <p — l. 
Recall that 

Skip) = E at 
1=1 

and that Sk{p) is the fcth elementary totally symmetric function in (71, ■ ■ ■ , g^(j,_iy 
By Proposition IT^ we infer that 

Sk{p) = Cp-i{k) (mod p). 

The latter congruence together with part 2 of Proposition El then yields the fol- 
lowing lemma. 
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Lemma 4 Let p be a prime. Then 



Skip) = ( T^-T-TT ) "^^^^-/l (mod p). 



This result is due to Moller ^3], who proved it in a different (and more round- 
about) way. Note that in case = 1 we have Si{p) = fi{p — 1) (mod p), a result 
first established by Gauss jH]. 



If ip{p — 1) < k, then Sk{p) = 0. If k = ip{p — 1), then we have, by Proposition 

Skip) = 9192 ■ ■■9^ip-i) = (-l)'^("-'^<fp-i(0) = l(mod p). 

If Lfip -1) > k, then Skip) = - 1) " ^) = (-l)^ap_i(A:)(mod p). 

Taking the results of these three cases together we infer that 

Skip) = (-l)^ap_i(A;)(mod p). 

The latter congruence and the congruence Skip) = Cp-iik) (mod p) are the start- 
ing point of our analysis of their value distribution which is carried out in the 
next two sections. 



9.2 Value distribution of Sk{p) mod p 

As we have seen Skip) = Cp-i(A;) (mod p) and so instead of the value distribution 
of Skip) ^od p we can investigate the value distribution of Cp_i(fc). Given a 
natural number k, let Dk = {qi, . . . ,qt} denote the set of prime divisors of k, 
ordered in size. We have t = i^Dk = uik) = J2p\k '^^^ following observation 
will play an important role. 

Proposition 20 On each prime p with vdSp ~ -'-) pi"escribed and fio^iP ~ 1) 
(respectively fio^iP " ^)'^) prescribed, Cp_i(/c) (respectively |cp_i(A;)|^ assumes the 
same value. If ^^kiP ~ 1) = 0? '^his value is zero. 

Our most general result regarding the value distribution of Cp_i(A;) reads as fol- 
lows. 

Theorem 5 Let t> 7^ 0. 

1) For every H > we have 

V 1 = c^ALiix) + O ( -4- ) , 

|cp_l(fc)l = f 

where the 0-constant depends at most on k and H and Cy is a rational number. 

2) // we assume in addition Conjecture 1, then the density of primes p with 
Cp_i(fc) = V exists and equals c^Ajl. 

Proof Part 1. Write i^OkiP ~ 1) = (^i, • • • ,eui{k))- If > ^qji^) + 2 for some 
1 < j < ^ik), then /i((p — l)/(p — 1, k)) = and hence Cp_i(A;) = 0. It follows 
there are only finitely many possibilities for vo^iP ~ 1) to be considered. To each 
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of them, by Proposition 1201 we can associate an unique non-zero value of Cp-i{k). 
We then have 

E 1= E -"^.(p-i)'' 

p<a; P'^x 
|cp_i(fc)|=t; ■'Dj.(P-l)eT 

for some effectively computable set T. The result then follows by Proposition IHl 
Part 2. The proof of part 2 is similar to that of part 1 and left to the reader. □ 

Conjecture 1 together with Theorem |21 (with r = 1 and k = 2) and the observation 
that Cp_i(l) = jJiip — 1), then yields the following result demonstrating part 2 of 
Theorem |31 (with k = l). 

Proposition 21 Assuming Conjecture 1 we have 

<5(cp_i(l)=j) = < I- A zfj = 0; 

Numerically we find up to the first 10^ primes 'densities' of 0.18732, 0.625881 and 
0.186799 respectively (for —1,0, 1). This should be compared to the conjectural 
values 0.1869 ■ ■ -, 0.6260 ■ ■ • and 0.1869 ■ ■ -, respectively. 

We demonstrate part 1 of Theorem El in Table 6 for k = 15, where we took x such 
that 7t{x) = 10^ and rounded the result in the column 'numerical' to the sixth 
decimal. 



Table 6: Density of values of |5'i5(p)| mod p 



V 


theoretical 


numerical 


approximate 





1 - 561A/475 


0.558339 


0.558178 


1 


9A/19 


0.177137 


0.177157 


2 


QA/19 


0.118091 


0.118138 


3 


2A/19 


0.039364 


0.039353 


4 


12A/95 


0.047237 


0.047316 


5 


12A/475 


0.009447 


0.009457 


8 


8A/95 


0.031491 


0.031508 


10 


8A/A75 


0.006298 


0.006317 


12 


8/1/285 


0.010497 


0.010486 


15 


8A/1425 


0.002100 


0.002090 



Let r be a prime and let k = r'^ for some e > 0. By Lemma 0] and (0) we deduce 

{/i(p-l) a Ur{p - I) = 0; 

r/(l-i)M^) ifz/,(p-l) = /, l</<e; (34) 
r>(^) ifi/,(p-l)>e. 

As regards to the average size of Cp-i{k) we have the following result. 
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Theorem 6 Let k be any natural number, and H any positive real number. Then 



q-1 



X 



\og^ X 



p<x q\k 

where the 0-constant depends at most on k and H . 

Corollary 3 The average of \cp-i{k)\/A over the primes p is a multiplicative 
function in k. 

Thus the average of |cp_i(/c)| is a semi-multiplicative function in k. The ana- 
log of this for the natural number case also holds true: Mn{\cn{m)\) is a semi- 
multiplicative function in m (we leave this as an exercise to the reader). 

Proof of Theorem^ Let us denote by S{q,e) the density of primes p such that 
z/q(p — 1) = e. By Lemma El one infers that 



6{q,e) 



1 - ^ if e = 0; 



if e > 1. 



By property 6 of Proposition El Proposition IHl the remark following Proposition 
inland we infer that the theorem holds with a constant given by 

6ius,ip - 1) ^ 0) n I 5(1, 0) + E ^(l^)5(l^ /) + l^'^'^^il, ^lik) + V 

q\k \ l<f<yq(k) 



On evaluating this constant, the result follows. 

Theorem mi is demonstrated in Table 7 (with again x such that tx{x) = 10^). 

Table 7: Average of |cp_i(A;)| 



□ 



k 


theoretical 


numerical 


approximate 


8 


AA 


1.494779 


1.495823 


21 


693A/205 


1.264572 


1.264153 


24 


36A/5 


2.689772 


2.692482 


27 


17A/5 


1.272214 


1.271450 


30 


126A/19 


2.479323 


2.479918 


36 


39A/5 


2.917172 


2.916855 



Our most general result in this section concerns the moments of Cp_i(/c). 
Theorem 7 Let k be a natural number. 

1). Let z 1 be a positive real number and H be a positive real number. Then 

p<x 
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where the 0-constant depends at most on k and H . 

2) . // we assume Conjecture 1 and j is an odd natural number, then we have 
X]p<x '^p-i(^)"' ~ <5(^(3^))j where the o-constant depends at most on k and j . 

Proof. Part 1. Proceeding as in the proof of Theorem IHl we infer that the result 
holds with constant 



q\k \ l<f<Vq{k) 

which is easily seen to equal the claimed constant. 

Part 2. Similar to the proof of Theorem El □ 

Remark. If in part 1 we take the limit for z tending to one, then we obtain 
Theorem El 

9.3 Value distribution of Sk{p) mod p 

Recall that Sk{p) = (— l)'^ap_i(/c)(mod p). It follows that \{sk{p) modp}| < 
\{an{k) : n > A presumably difficult question is to investigate under which 
conditions equality holds here. A related question is what the set {ap-i(/c) : k > 
1, p prime } looks like (cf. Suzuki's result in Section 4.1). 

First we consider the value distribution of Sk{p) for some small values of k. When 
= 1 we have Si{p) = Si{p) and hence for this case we refer to the previous 
section. Note that S2{p) = Ylii<i<j<^p{p-i)9i9j- remainder of this section 

we assume that p > 2 (some of the assertions will also be valid for p = 2). Using 
Lemma m we infer that 

S2{P) = liS.ipf - S,ip)) ^ \ (^(^(p - 1))^ - ^^-"(^)) (mod p)- 

where the latter identity is Theorem 3 of 0. By (0) we have 

(p{p — 1) _ / 2 iip= l(mod 4); 



ip{^) \1 ifp = 3(mod4). 

On using this we infer that 

= ifp^l(mod4); , . 

U(p-l)(/i(p-l) + l)/2 ifp = 3(mod4). ^ ' 

(For notational convenience we will write Sk{p) = a for Sk{p) = a(mod p).) Since 
82(2) = it follows that the reduction of S2{p) mod p is in {—1,0, 1}. The fol- 
lowing assertion conjecturally resolves the question of the Dence brothers stated 
in the introduction. 
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Proposition 22 We have, assuming the validity of Conjecture 1, 

( 4 ^fj = -1/ 

6{s2{p)^j) = { I- A ifj = 0; 

[ if J = 1. 

Proof. Equation (jH^j) suggests to consider the primes p with p = l(mod 4) and 
p = 3 (mod 4) separately. Note that for x > 2, 

p<x. p=3(mod 4) p<x 

Ai(p-l) = l H(P-1) = 1 

The latter quantity we (conditionally) evaluated in Proposition |^ It follows 
that the density of primes p such that p = 3 (mod 4) and S2{p) = 1 equals A/2 
and that the density of primes p such that p = 3(mod 4) and S2(p) = equals 
1/2 - A/2. 

The case p = l(mod 4) can be dealt with Proposition IHl (with t = 1, qi = 2 
and ei = 2) and Conjecture 1. The results are summed up in Table 8. □ 



Table 8: Value distribution of S2{p) mod p 



5 


-1 





1 


+ 


U2{p - 1) < 1 

numerical 

7r(x) = 10^ 




0.000000 
0.000000 


1/2 - A/2 
0.313022 
0.313403 


A/2 
0.186978 
0.186798 


1/2 
0.500000 
0.500201 


i^2ip-l) = 2 
numerical 

tt{x) = 10^ 


A/ A 
0.093489 
0.093939 


1/2 - A/2 
0.313022 
0.312813 


A/ A 
0.093489 
0.093047 


1/2 
0.500000 
0.499799 


+ 

numerical 

tt{x) = 10^ 


A/ A 
0.093489 
0.093939 


1-A 
0.626044 
0.626216 


3A/A 
0.280467 
0.279845 


1 

1.000000 
1.000000 



By Newton's formula Sk{p) can be expressed as a polynomial in Sr{p) with 1 < 
r < k. To be precise, we have 

where the sum extends over all solutions {ki, k2, ■ ■ ■) of ki + 2k2 + ■ ■ ■ = k. Next 
we consider s^^p). Taking = 3 in (jHUI) we compute 

, , S,{pr + 2S,{p)-3S^{p)S2{p) 

Invoking Lemma |3 a more explicit formula for ss{p) can then be derived. With 

f3 = u^ip — 1) we find 

rMp-l)(Mp-l) + l)/2 if/? = 0; 

S3(p)= < Mp-i)(Mp-i)-i)/2 if/? = i; 

In each of these cases Proposition IHl and Conjecture 1 yield the densities with 
which the values —1,0 and 1 are assumed. The results are summed up in Table 
9. 
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Table 9: Value distribution of s^i^p) 






1 

— i 


U 


i 




/3 = 
numerical 

7r[X) = iU 




0.000000 

U.UUUUUU 


1/2 - 3A/10 
0.387813 
U.oooioo 


3A/10 
0.112187 

U.iizUoO 


1/2 
0.500000 

U.oUUi 1 


(3 = 1 
numerical 

TTyX ) — iU 




0.000000 
n nnnnnn 

U.UUUUUU 


1/3 - A/5 
0.258542 

U./Ooo ( o 


A/5 

0.074791 
n n7/i ssQ 

U.U ( 4ooo 


1/3 
0.333333 

U.OOOZOi 


/3 > 2 
numerical 
tt{x) = 10^ 


A/15 
0.024930 
0.025031 


1/6 - 2A/15 
0.116806 
0.116729 


A/15 
0.024930 
0.024809 


1/6 
0.166666 
0.166569 


+ 

numerical 

7t{x) = 10^ 


A/15 
0.024930 
0.025030 


1 - 19A/30 
0.763161 
0.763243 


17A/30 
0.211908 
0.211727 


1 

1.000000 
1.000000 



Adding up the various contributions we infer the following result. 
Proposition 23 We have, assuming the validity of Conjecture 1, 

{A if j = _i- 

i - §A if 3 = 0; ' 

Similarly, for S4(p) we find, with a = h'2{p — 1) and P = u^i^p — 1), 

S4{p) = 



n{p-l){li{p~l) + l)/2 if a = 1 and /5 = 0; 
fj,{p-l)ll- fi{p-l))/2 if a = 1 and/3 > 1; 

M¥)(M¥) + l)/2 if« = 2; 



Proceeding as before we infer the following result. 

Proposition 24 We have, assuming the validity of Conjecture 1, 

r ^fj = -1; 

6{s,{p)^3) = { I- A if 3 = 0; 

I ^/^- = 1- 

The above approach of dealing with the value distribution of Sfc(p)(modp) for 
small p clearly becomes more laborious as k increases. The following result is 
more systematic; it is analogous to Theorem 0] 

Theorem 8 Let k > 1 be fixed. Let Mk = k np</c P ^'^^ let r > k be any prime. 
Then 

Ax T-r q{q — 2) sr^ ad(k) + adr{k) Yj q — 1 x 



E(k] — TT ~ ^ ^dyi^j -r Udryi^) -TT 

"^"^^ ^ ~ 21ogx J-J- g2 1 2^ ^ llg_2 ' "^logx' 

p<x 2<q<k ^ ^ d\M^ q\d ^ ^ 

2\d q>2 
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Ifvy^O, then 



log X loe X ' 

tlp_l(fc)=u 



where 



«^(*) 4 n #^ 



2<g<fc ^ ^ 



/ 1 _i 1-1^ 

d|A/^, 2\d q\d ^ ti|Mj., 2|d q\d ^ ; 

\ «d('=)=i' 1>2 a^^{k)=v q>2 / 



In both cases the o-constant depends at most on k. 

Proof. Put Nk = lcm(l, 2, - ■ ■ ,k) Y[p<kP- ^et S{k) = {qi, . . . , ^^^(fc)} be the set 
of primes not exceeding k and r2 = Nf^ Y[p<kP- Furthermore, let d\Nk. We first 
consider ^p<^ (p-i r2)=d/^((^' ~ ^)/^)- Let q G S{k) and suppose that p satisfies 
{p — 1, = (i. Since i^q{r2) > i^q{d), it follows that Ug^p — 1) = Pg^d). Note that 
fj,{{p — l)/d) = fis(k){p — !)• Conjecture 1 implies now that 

" ^ /^5(fc)(p- 1) = (37) 

p<x, {p-l,r2)=d P<^ 

where the o-constant depends at most on k. 
Next we consider 

E '"^^)'= E (38) 

p<a; p<a; 
(p-l,r2) = ti '^S(fe)(P~l)^(''9l('^)'---'''9^(fc) 

By Proposition El we then find that 



^ ^ ^„<7V 9-1/ g(g-l)/loga; 'logo;' 



p<x q<k \ J / q^^j 

(p— l,r) — d q\d 



Using Proposition 13 we have 

J2a,^,{k) = j2 E UiidM^r + BMK^)) ■ 

p<x d\Nk P<^ ^ ^ 

(p-l.'-2)=d 

From the latter formula in combination with (j37|) and (j38|) . we infer that 

p<x ^ 2<g<fc ^ ^ d|iVj. g|d ^ ^ 

2|d >J>2 

On noting that Ai{d) = {ad{k) + adr{k))/2 and using (jTU)), the proof of the first 
assertion is completed. The proof of the second assertion is a variation of the 
argument above and left to the reader. □ 
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The latter theorem is demonstrated in Table 10. Using that the sum of the den- 
sities equals one, one can easily find the density of primes for which ap_i(fc) = 
from the table, so we omitted this information. Notice that the averages given for 
1 < k < 4 are the same as can be inferred from the earlier presented results on 
the value distribution of Si{p) , S2{p) , s^i^p) and s^i^p) respectively. The averages 
given also seem to be consistent with numerical simulations. 



Table 10: Conjectural value distribution of 6{ap_i{k) = v)/A 



J 


-2 


-1 


1 


2 


Average 


k = 1 





1/2 


1/2 








k = 2 





1/4 


3/4 





1/2 


k = 3 





17/30 


1/15 





-1/2 


k = 4 





13/40 


27/40 





7/20 


k = 5 





69/90 


6/95 





-3/10 


k = 6 





443/1140 


47/95 





121/1140 


k = 7 





13989/54530 


358/1435 


24/3895 


1/190 


k = 8 





16703/62320 


35873/62320 





1917/6232 


k = 9 





31477/70110 


2129/35055 





-9073/23370 


k=10 





267/820 


505/1558 





-23/15580 



It is not difficult to see that 6{ap_i{k) = f ) > if and only if ap-i{k) = v for some 
prime p. We infer from the table that there is no prime p such that ap_i(7) = —2, 
whereas by Table 4 it follows that (5(a„(7) = —2) = 0.001055 ■■ ■. Indeed, the 
example following Proposition ^2 shows that there are no even integers n for 
which a„(7) = —2. 



10 Open problems 

During research usually more questions are being raised than being resolved and 
the present work is no exception. We conclude by formulating some questions. 
Gallot's data show that 2, 5 and 9 hold true for k < 100. 

1) Is Conjecture 1 true ? 

2) Is MoUer's conjecture that < < 1/2 true ? 

3) What is the behaviour of ek for large k ? 

4) Compute M(efc). 

5) Is efc/cnp<,b + l)eZ? 

6) Determine {ap_i(fc) : p prime, k G N}. 

7) Does an{k) = —v, v < 0, imply there is an integer m with am{k) = v 1 

8) Is \B{k)\BQ{k)\ > 1 infinitely often ? 

9) Is 5{aJ{k) = 1) > 5{anik) = -1) ? 

10) What is the behaviour of 6{an{k) ^ {—1, 0, 1}) for large k ? 

Acknowledgement. We are very grateful to dr. Gallot for assistance in com- 
puting some of the tables. We thank Prof. T. Dence for sending [S] and Prof. 
Thangadurai for sending us his beautiful survey paper Prof. Bachman for 
some helpful e-mail correspondence and dr. Tegelaar for making some plots using 
Mathematica. 
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Table 11: Value distribution of a„(/c) for 1 < A; < 61 

(The number between brackets is kY[p<:k{p + ^)^k-) 

ei - 0(0.000000) [0] 
Vi[-1] = 1/2(0.500000) 
Vi[+1] = 1/2(0.500000) 



62 
V2 
V2 

es 
Vs 

64 
^4 
^4 

65 

66 
^6 
^6 

67 

Vr 

Vr 
V7 
V7 

68 

^8 

Vs 

69 
^9 



= 1/2(0.500000) [3] 
-1] - 1/12(0.083333) 
+1] = 7/12(0.583333) 

= 1/6(0. 166667)[6] 
-1] = 5/24(0.208333) 
+1] = 3/8(0.375000) 

= 1/3(0.333333) [16] 
-1] = 1/6(0.166667) 
+1] = 1/2(0.500000) 

= 1/8(0.125000) [45] 

1] = 13/80(0.162500) 
+1] = 23/80(0.287500) 

= 7/24(0.291667) [126] 

1] = 25/144(0.173611) 
+1] = 67/144(0.465278) 

= 1/18(0.055556) [224] 
-2] = 1/576(0.001736) 

-1] = 577/2688(0.214658) 
+1] = 731/2688(0.271949) 
+2] = 1/1152(0.000868) 

= 7/24(0.291667) [1344] 

-1] = 1/8(0.125000) 
+1] = 5/12(0.416667) 

= 19/144(0. 131944) [684] 
-1] = 65/384(0.169271) 
+1] = 347/1152(0.301215) 



do = 31/160(0.193750)[1116] 
Vio[-l] = 161/960(0.167708) 
Viol+1] = 347/960(0.361458) 
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= 1/16(0.062500) [4752] 
-2] = 1/2304(0.000434) 
-1] = 8299/50688(0.163727) 
+1] = 11489/50688(0.226661) 
+2] = 1/4608(0.000217) 

= 55/192(0.286458) [23760] 
-1] = 349/2304(0.151476) 
+ 1] = 1009/2304(0.437934) 

= 13/288(0.045139) [56784] 
-2] = 43/48384(0.000889) 
-1] = 219269/1257984(0.174302) 
+ 1] = 277171/1257984(0.220330) 
+2] = 43/96768(0.000444) 

= 61/288(0.211806) [286944] 
-2] = 13/21504(0.000605) 
-1] = 2395/21504(0.111375) 
+ 1] = 2319/7168(0.323521) 
+2] = 1/2304(0.000434) 

= 2287/20160(0.113442) [164664] 
-2] = 13/32256(0.000403) 
-1] = 1345/7168(0.187640) 
+ 1] = 97247/322560(0.301485) 
+2] = 13/64512(0.000202) 

= 733/4032(0.181796) [281472] 
-2] = 5/21504(0.000233) 

-1] = 12149/64512(0.188322) 
+ 1] = 1127/3072(0.366862) 
+2] = 5/2688(0.001860) 

= 667/8064(0.082713) [2449224] 

-3] = 5/580608(0.000009) 

-2] = 281/193536(0.001452) 

-1] = 2353487/19740672(0.119220) 

+1] = 1981753/9870336(0.200779) 

+2] = 197/96768(0.002036) 

+3] = 5/1161216(0.000004) 
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= 79/336(0.235119) [7371648] 
-2] = 961/1161216(0.000828) 
-1] = 5575/43008(0.129627) 
+1] = 212369/580608(0.365770) 
+2] = 341/1161216(0.000294) 
+3] = 17/1161216(0.000015) 

= 55/1344(0.040923) [27086400] 
-3] = 19/2322432(0.000008) 
-2] = 67813/69672960(0.000973) 
-1] = 42731243/264757248(0.161398) 
+ 1] = 67086449/330946560(0.202711) 
+2] = 54641/69672960(0.000784) 
+3] = 19/4644864(0.000004) 

= 221/960(0.230208) [160392960] 
-3] = 1/221184(0.000005) 
-2] = 307/165888(0.001851) 
-1] = 1417037/11612160(0.122030) 
+ 1] = 8240789/23224320(0.354834) 
+2] = 361/645120(0.000560) 

= 8207/120960(0.067849) [49635936] 
-3] = 1/331776(0.000003) 
-2] = 138259/69672960(0.001984) 
-1] = 14659501/69672960(0.210404) 
+ 1] = 69503/248832(0.279317) 
+2] = 101363/69672960(0.001455) 
+3] = 1/663552(0.000002) 

= 8467/95040(0.089089) [68277888] 

-3] = 1/221184(0.000005) 

-2] = 7/18432(0.000380) 

-1] = 8101001/42577920(0.190263) 

+1] = 70728809/255467520(0.276860) 

+2] = 3751/2322432(0.001615) 

+3] = 19/1658880(0.000011) 

= 629/11520(0.054601) [1049956992] 
-4] = 1/7962624(0.000000) 

-3] = 107/7962624(0.000013) 

-2] = 578371/418037760(0.001384) 

-1] = 5016105893/38459473920(0.130426) 

+1] = 1789013287/9614868480(0.186067) 

+2] = 730129/836075520(0.000873) 

+3] = 107/15925248(0.000007) 

+4] = 1/15925248(0.000000) 
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624 = 7327/24192(0.302869) [6077306880 
-3] = 853/69672960(0.000012) 
-2] = 247859/278691840(0.000889) 
-1] = 898117/9289728(0.096679) 
+ 1] = 55697567/139345920(0.399707) 
+2] = 8227/10321920(0.000797) 
+3] = 713/34836480(0.000020) 



^24 
^24 
^24 
V2, 
V2, 
V2^ 

625 
^25 
1^25 
V^25 
V25 
F25 
1^25 

626 
^26 
1^26 
^26 
V^26 
V^26 
^26 

627 
V^27 
V^27 
^27 
1^27 
V27 
V27 

628 
1^28 
^28 
V2S 
^28 
^28 
T^28 
T^28 

629 
V29 
V2^ 

V29 
^29 



= 1087/18144(0.059910) [1252224000] 

-3] = 25/7962624(0.000003) 

-2] = 81607/238878720(0.000342) 

-1] = 1375923029/8360755200(0.164569) 

+ 1] = 469535011/2090188800(0.224638) 

+2] = 88489/334430208(0.000265) 

+3] = 25/15925248(0.000002) 

= 234433/1572480(0. 149085) [3240801792] 

-3] = 25/5308416(0.000005) 

-2] = 11149/15482880(0.000720) 

-1] = 837724879/7245987840(0.115612) 

+ 1] = 127807265/483065856(0.264575) 

+2] = 54667/69672960(0.000785) 

+3] = 1289/557383680(0.000002) 

= 33491/362880(0.092292) [2083408128] 
-3] = 811/185794560(0.000004) 
-2] = 311893/836075520(0.000373) 
-1] = 244503443/1672151040(0.146221) 
+1] = 399281329/1672151040(0.238783) 
+2] = 201713/836075520(0.000241) 
+3] = 139/61931520(0.000002) 

= 84047/483840(0.173708) [4066530048] 
-3] = 1073/69672960(0.000015) 

-2] = 202387/278691840(0.000726) 
-1] = 5770579/39813120(0.144942) 
+ 1] = 88832659/278691840(0.318749) 
+2] = 96629/139345920(0.000693) 
+3] = 1/241920(0.000004) 
+4] = 1/7741440(0.000000) 

= 5021/103680(0.048428) [35225729280] 
-4] = 1837/25082265600(0.000000) 
-3] = 215897/16721510400(0.000013) 
-2] = 56528869/50164531200(0.001127) 
-1] = 5288720741/41564897280(0.127240) 



46 



V29 
V29 
^29 

eso 

Vm 
V^o 
Vm 
Vso 
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V^31 
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V^31 
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^33 

^33 
V^33 
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+1] = 256004923169/1454771404800(0.175976) 
+2] = 48779039/50164531200(0.000972) 
+3] = 24511/1857945600(0.000013) 
+4] = 3883/50164531200(0.000000) 

= 45893/241920(0. 189703) [142745587200] 
-3] = 55183/8360755200(0.000007) 
-2] = 2194267/2090188800(0.001050) 
-1] = 541742161/3344302080(0.161990) 
+ 1] = 5867204267/16721510400(0.350878) 
+2] = 34039/23887872(0.001425) 
+3] = 450059/16721510400(0.000027) 
+4] = 8107/8360755200(0.000001) 
+5] = 29/5573836800(0.000000) 

= 155/5376(0.028832) [717382656000] 

-4] = 341/7644119040(0.000000) 

-3] = 25652623/802632499200(0.000032) 

-2] = 1728495157/1605264998400(0.001077) 

-1] = 4263256278479/24881607475200(0.171342) 

+ 1] = 4956093793291/24881607475200(0.199187) 

+2] = 2463624983/1605264998400(0.001535) 

+3] = 44559467/802632499200(0.000056) 

+4] = 341/15288238080(0.000000) 

= 4381/17920(0.244475) [6279166033920] 
-4] = 143/5096079360(0.000000) 
-3] = 143131/19818086400(0.000007) 
-2] = 32276503/38220595200(0.000844) 
-1] = 45508919251/535088332800(0.085049) 
+ 1] = 24918610681/76441190400(0.325984) 
+2] = 667940893/267544166400(0.002497) 
+3] = 8249729/107017666560(0.000077) 
+4] = 1207481/178362777600(0.000007) 

= 294509/3193344(0.092226) [2442775449600] 
-4] = 20065/10701766656(0.000002) 
-3] = 80453581/802632499200(0.000100) 
-2] = 412044869/401316249600(0.001027) 
-1] = 2239952947681/17657914982400(0.126853) 
+1] = 240607474937/1103619686400(0.218017) 
+2] = 1136159459/802632499200(0.001416) 
+3] = 309052627/1605264998400(0.000193) 
+4] = 398297/107017666560(0.000004) 
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634 = 268801/3525120(0.076253) [2080906813440] 



1^34 
^34 

^34 
V^34 
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-4] = 143/5096079360(0.000000) 

-3] = 20977/5096079360(0.000004) 

-2] = 1313801/1189085184(0.001105) 

-1] = 1500096902327/9096501657600(0.164909) 

+ 1] = 135360848387/568531353600(0.238089) 

+2] = 80945987/33443020800(0.002420) 

+3] = 3130531/21403533312(0.000146) 

+4] = 1071667/267544166400(0.000004) 

+5] = 719/19818086400(0.000000) 

= 69809/870912(0.080156) [2251759104000] 
-5] = 323/2090188800(0.000000) 
-4] = 242639/25082265600(0.000010) 
-3] = 280214099/1605264998400(0.000175) 
-2] = 428425901/229323571200(0.001868) 
-1] = 15515446957/107017666560(0.144980) 
+ 1] = 120969894043/535088332800(0.226075) 
+2] = 2211435001/1605264998400(0.001378) 
+3] = 292252153/1605264998400(0.000182) 
+4] = 4691237/321052999680(0.000015) 
+5] = 323/1045094400(0.000000) 

= 41333/207360(0.199330) [5759584911360] 
-3] = 5072809/76441190400(0.000066) 

-2] = 31579379/26754416640(0.001180) 
-1] = 92196104039/535088332800(0.172301) 
+ 1] = 12419782037/33443020800(0.371371) 
+2] = 124255577/107017666560(0.001161) 
+3] = 78450899/535088332800(0.000147) 
+4] = 22711/1651507200(0.000014) 
+5] = 56381/178362777600(0.000000) 
+6] = 14633/535088332800(0.000000) 

= 19073/544320(0.035040) [39542741729280] 

-5] = 30631/3210529996800(0.000000) 

-4] = 56904067/61000069939200(0.000001) 

-3] = 60188399/1742859141120(0.000035) 

-2] = 76311280099/61000069939200(0.001251) 

-1] = 323928333998587/2257002587750400(0.143521) 

+1] = 80619151138297/451400517550080(0.178598) 

+2] = 73686178829/61000069939200(0.001208) 

+3] = 763875323/15250017484800(0.000050) 

+4] = 51705001/30500034969600(0.000002) 

+5] = 30631/1605264998400(0.000000) 
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= 246539/1451520(0.169849) [196855040655360] 
-4] = 7955363/2541669580800(0.000003) 
-3] = 18756869/150617456640(0.000125) 
-2] = 1660704037/1694446387200(0.000980) 
-1] = 137023551787/1694446387200(0.080866) 
+ 1] = 1265500142927/5083339161600(0.248951) 
+2] = 4927088377/2541669580800(0.001939) 
+3] = 96190153/1270834790400(0.000076) 
+4] = 17066753/10166678323200(0.000002) 
+5] = 629/356725555200(0.000000) 

= 961879/12804480(0.075121) [89355942789120] 

-6] = 1/179159040(0.000000) 

-5] = 179/627056640(0.000000) 

-4] = 77322241/12200013987840(0.000006) 

-3] = 6774345563/61000069939200(0.000111) 

-2] = 2729296757/2033335664640(0.001342) 

-1] = 121853678359577/793000909209600(0.153661) 

+1] = 181373355248401/793000909209600(0.228718) 

+2] - 13458490453/10166678323200(0.001324) 

+3] = 1683824087/12200013987840(0.000138) 

+4] = 7100759/642105999360(0.000011) 

+5] = 179/313528320(0.000001) 

+6] = 1/89579520(0.000000) 

= 3360923/18385920(0.182799) [223014717849600] 
-5] = 198679/26754416640(0.000007) 
-4] = 2961557/4066671329280(0.000001) 
-3] = 5036034173/20333356646400(0.000248) 
-2] = 30381984071/20333356646400(0.001494) 
-1] = 2957102734199/20333356646400(0.145431) 
+1] = 3360988664941/10166678323200(0.330589) 
+2] = 6478070507/10166678323200(0.000637) 
+3] = 292144933/6777785548800(0.000043) 
+4] = 9938893/5083339161600(0.000002) 
+5] = 4597979/20333356646400(0.000000) 
+6] = 13/30576476160(0.000000) 

= 7313311/165473280(0.044196) [2321237890621440] 
-6] = 61468511/512400587489280(0.000000) 

-5] = 905963641/183000209817600(0.000005) 

-4] = 739141561/128100146872320(0.000006) 

-3] = 95451411817/427000489574400(0.000224) 

-2] = 184063646579/73200083927040(0.002515) 

-1] = 3161605893907487/26260530108825600(0.120394) 

+1] = 8685630577953149/52521060217651200(0.165374) 

+2] = 2865735759811/1281001468723200(0.002237) 
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+3] = 65007805063/427000489574400(0.000152) 
+4] = 12935415793/2562002937446400(0.000005) 
+5] = 46277477/18300020981760(0.000003) 
+6] = 77476183/1281001468723200(0.000000) 

= 6123037/27578880(0.222019) [11945083135426560] 
-6] = 396751/30500034969600(0.000000) 
-5] = 14830471/20333356646400(0.000001) 
-4] = 1896873157/142333496524800(0.000013) 
-3] = 3201729511/13343765299200(0.000240) 
-2] = 33789705539/18977799536640(0.001780) 
-1] = 93023983380571/854000979148800(0.108927) 
+1] = 283373054708011/854000979148800(0.331818) 
+2] = 705994382021/427000489574400(0.001653) 
+3] = 6305790613/122000139878400(0.000052) 
+4] = 338310517/427000489574400(0.000001) 
+5] = 159624961/854000979148800(0.000000) 
+6] = 7585/1265186635776(0.000000) 

= 1241881 /47278080(0.026268) [63663558680494080] 

-7] = 899/20600900812800(0.000000) 

-6] = 239374189/56364064623820800(0.000000) 

-5] = 8097232133/56364064623820800(0.000000) 

-4] = 35174564923/12525347694182400(0.000003) 

-3] = 38747727775/501013907767296(0.000077) 

-2] = 34634252980079/18788021541273600(0.001843) 

-1] = 779722198548416059/4847309557648588800(0.160857) 

+1] = 1906106119849907/10204862226628608(0.186784) 

+2] = 8259199022633/4175115898060800(0.001978) 

+3] = 536250514949/5368006154649600(0.000100) 

+4] = 2687713733/751520861650944(0.000004) 

+5] = 6777068293/56364064623820800(0.000000) 

+6] = 281711701/112728129247641600(0.000000) 

+7] = 899/41201801625600(0.000000) 

= 27117469/151683840(0.178776) [443368483462840320] 
-7] = 19499/1569592442880(0.000000) 
-6] = 483086393/521889487257600(0.000001) 
-5] = 11256275861/596445128294400(0.000019) 
-4] = 664972869923/7515208616509440(0.000088) 
-3] = 1013158170881/12525347694182400(0.000081) 
-2] = 35958299894189/12525347694182400(0.002871) 
-1] = 209668413473641/2087557949030400(0.100437) 
+1] = 279518297968507/988843239014400(0.282672) 
+2] = 2998032785137/2087557949030400(0.001436) 
+3] = 1257623773787/37576043082547200(0.000033) 
+4] = 64724291587/37576043082547200(0.000002) 
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+5] = 50451869/1977686478028800(0.000000) 
+6] = 65231/357867076976640(0.000000) 

= 3979615/45505152(0.087454) [221817245368320000] 

-7] = 33263/18788021541273600(0.000000) 
-6] = 566438617/10248011749785600(0.000000) 
-5] = 269113485517/112728129247641600(0.000002) 
-4] = 12278820607/4697005385318400(0.000003) 
-3] = 151929636707/1503041723301888(0.000101) 
-2] = 46254490783673/28182032311910400(0.001641) 
-1] = 6148678815334457/37576043082547200(0.163633) 
+1] = 295732896088667/1174251346329600(0.251848) 
+2] = 149807819591587/112728129247641600(0.001329) 
+3] = 2209747042441/37576043082547200(0.000059) 
+4] = 340486901/191714505523200(0.000002) 
+5] = 778779179/644160738557952(0.000001) 
+6] = 314248783/11272812924764160(0.000000) 
+7] = 33263/37576043082547200(0.000000) 

= 1634903647/27909826560(0.058578) [151878101868380160] 

-7] = 593147/1404606873600(0.000000) 

-6] = 8041/200658124800(0.000000) 

-5] = 461160835213/37576043082547200(0.000012) 

-4] = 87404403091/18788021541273600(0.000005) 

-3] = 772557276649/18788021541273600(0.000041) 

-2] = 14774599478917/6262673847091200(0.002359) 

-1] = 30866961312674209/172849798179717120(0.178577) 

+ 1] = 1461674766295439/6173207077847040(0.236777) 

+2] = 22857774451559/9394010770636800(0.002433) 

+3] = 533933483887/4175115898060800(0.000128) 

+4] = 218641508917/18788021541273600(0.000012) 

+5] = 2284159373/2087557949030400(0.000001) 

+6] = 147831841/2684003077324800(0.000000) 

+7] = 8611853/7515208616509440(0.000000) 

+8] = 19/2809213747200(0.000000) 

= 64767749/1213470720(0.053374) [6786891421192028160] 

-9] = 1147/715734153953280(0.000000) 

-8] = 56233703/169092193871462400(0.000000) 

-7] = 8196769949/676368775485849600(0.000000) 

-6] = 231016391/27606888795340800(0.000000) 

-5] = 115585824097/216438008155471872(0.000001) 

-4] = 32560081364423/2705475101943398400(0.000012) 

-3] = 31938251933447/193248221567385600(0.000165) 

-2] = 10645199114119417/5410950203886796800(0.001967) 

-1] = 27034488309751213037/254314659582679449600(0.106303) 

+1] = 115130557033516291/722484828359884800(0.159354) 
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+2] = 5876022667634491/2705475101943398400(0.002172) 
+3] = 31091772610243/216438008155471872(0.000144) 
+4] = 3905889918467/541095020388679680(0.000007) 
+5] = 191742884821/676368775485849600(0.000000) 
+6] = 3669340891/772992886269542400(0.000000) 
+7] = 32823296417/5410950203886796800(0.000000) 
+8] = 56233703/338184387742924800(0.000000) 
+9] = 1147/1431468307906560(0.000000) 

= 3288940867/12134707200(0.271036) [35197477704016330752] 

-7] = 25505321/94928950945382400(0.000000) 

-6] = 8172767/645775176499200(0.000000) 

-5] = 1007637241/47464475472691200(0.000000) 

-4] = 194589824999/163968187996569600(0.000001) 

-3] = 23433243267407/901825033981132800(0.000026) 

-2] = 181000587404083/150304172330188800(0.001204) 

-1] = 174759438479524523/1803650067962265600(0.096892) 

+1] = 109207543272255949/300608344660377600(0.363288) 

+2] = 5722102669224553/1803650067962265600(0.003173) 

+3] = 42223142910659/180365006796226560(0.000234) 

+4] = 4708918459513/257664295423180800(0.000018) 

+5] = 18463193263/9394010770636800(0.000002) 

+6] = 23339473/303644792586240(0.000000) 

+7] = 1193860483/66801854368972800(0.000000) 

+8] = 8671/10844456878080(0.000000) 

+9] = 519961/54656062665523200(0.000000) 

= 10989331/260029440(0.042262)[5602583726551203840] 

-7] = 131153887/901825033981132800(0.000000) 

-6] = 8877334111/676368775485849600(0.000000) 

-5] = 321058535947/676368775485849600(0.000000) 

-4] = 188187460081/17568020142489600(0.000011) 

-3] = 149860268269/1294485694709760(0.000116) 

-2] = 9555417316875511/5410950203886796800(0.001766) 

-1] = 77997264836317111/541095020388679680(0.144147) 

+1] = 12026849533501/64724284735488(0.185817) 

+2] = 5202805469811217/2705475101943398400(0.001923) 

+3] = 106923208106387/541095020388679680(0.000198) 

+4] = 24986800232473/1352737550971699200(0.000018) 

+5] = 2087822266097/2705475101943398400(0.000001) 

+6] = 87856669279/5410950203886796800(0.000000) 

+7] = 8813377/64416073855795200(0.000000) 

= 379369/2996224(0.126616) [17127781233131520000] 
-7] = 2101783/47464475472691200(0.000000) 

-6] = 164943659/33400927184486400(0.000000) 

-5] = 465186406151/1803650067962265600(0.000000) 
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-4] = 49542317623/9109343777587200(0.000005) 
-3] = 6825915940891/100202781553459200(0.000068) 
-2] = 3811750874427449/1803650067962265600(0.002113) 
-1] = 29887357972193783/225456258495283200(0.132564) 
+1] = 6548160578576831/25050695388364800(0.261396) 
+2] = 708473647746719/901825033981132800(0.000786) 
+3] = 356805154396393/1803650067962265600(0.000198) 
+4] = 11564370480493/901825033981132800(0.000013) 
+5] = 893786131687/225456258495283200(0.000004) 
+6] = 468039803071/1803650067962265600(0.000000) 
+7] = 104062907/25050695388364800(0.000000) 
+8] = 59551/12269728353484800(0.000000) 

= 12998582579/206290022400(0.063011) [8694237351200292864] 
-8] = 302799841/71196713209036800(0.000000) 

-7] = 200198636767/1082190040777359360(0.000000) 

-6] = 4232845423/28478685283614720(0.000000) 

-5] = 296867080061/110427555181363200(0.000003) 

-4] = 28279018507781/2705475101943398400(0.000010) 

-3] - 5285026941479/48312055391846400(0.000109) 

-2] = 2400296876091767/1352737550971699200(0.001774) 

-1] = 36173757934527517/238925073937858560(0.151402) 

+1] = 9847006606815082297/45993076733037772800(0.214098) 

+2] = 9792585533485231/5410950203886796800(0.001810) 

+3] = 67216527293609/386496443134771200(0.000174) 

+4] = 105773214248153/5410950203886796800(0.000020) 

+5] = 1496553336277/284786852836147200(0.000005) 

+6] = 1592766194951/5410950203886796800(0.000000) 

+7] = 10213095919/27606888795340800(0.000000) 

+8] = 302799841/35598356604518400(0.000000) 

= 4136435441/39437798400(0.104885) [14755722292396818432] 

-7] = 154583/1078738078924800(0.000000) 

-6] = 8275709/89894839910400(0.000000) 

-5] = 5082975040331/601216689320755200(0.000008) 

-4] = 4086365658149/200405563106918400(0.000020) 

-3] = 146925185217599/450912516990566400(0.000326) 

-2] = 3225891061904431/901825033981132800(0.003577) 

-1] = 3629547434141844893/23447450883509452800(0.154795) 

+1] = 6208463813049931387/23447450883509452800(0.264782) 

+2] = 2533165683241/2109532243230720(0.001201) 

+3] = 188653240419367/901825033981132800(0.000209) 

+4] = 3227651747927/112728129247641600(0.000029) 

+5] = 3620914111/1937325529497600(0.000002) 

+6] = 233961671281/1803650067962265600(0.000000) 

+7] = 3363306073/150304172330188800(0.000000) 

+8] = 35069/28757175828480(0.000000) 
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V^2 
^52 

653 
^53 
^53 
^53 
V^53 
^53 
V^53 
V^53 
^53 
V^53 
V53 
^53 
^53 
^53 
V^53 
^53 
^53 
^53 
^53 
V^53 
^53 
^53 
V^53 
V^53 
^53 
^53 
T^53 

654 

V,, 

1^54 

1^54 
1^54 
V,, 
^54 
T^54 
1^54 
V54 
V,, 
^54 
^54 
1^54 
V^54 



+9] = 218671/7910745912115200(0.000000) 
+10] = 10013/47464475472691200(0.000000) 

= 1753484317/36404121600(0.048167) [372961927410354487296] 

-13] = 1/4831205539184640(0.000000) 

-12] = 137923/3043659489686323200(0.000000) 

-11] = 70055/15218297448431616(0.000000) 

-10] = 24043013971/97397103669962342400(0.000000) 

-9] = 536992451/81984093998284800(0.000000) 

-8] = 495270318673/7304782775247175680(0.000000) 

-7] = 960637729199/24349275917490585600(0.000000) 

-6] = 214231483919113/292191311009887027200(0.000001) 

-5] = 56397449178799/13281423227722137600(0.000004) 

-4] = 3061411267633/216438008155471872(0.000014) 

-3] = 15076596535394231/97397103669962342400(0.000155) 

-2] = 250559460613448869/146095655504943513600(0.001715) 

-1] = 212226881850682202239/1935767435440501555200(0.109634) 

+1] = 243813332843085582721/1548613948352401244160(0.157440) 

+2] = 105517436829081911/58438262201977405440(0.001806) 

+3] = 9720451206859513/48698551834981171200(0.000200) 

+4] = 332885526845159/16232850611660390400(0.000021) 

+5] = 2199855412954037/292191311009887027200(0.000008) 

+6] = 393376392231287/292191311009887027200(0.000001) 

+7] = 1765246803409/24349275917490585600(0.000000) 

+8] = 10420346053/76892450265759744(0.000000) 

+9] = 536992451/40992046999142400(0.000000) 

+10] = 24043013971/48698551834981171200(0.000000) 

+11] = 70055/7609148724215808(0.000000) 

+12] = 137923/1521829744843161600(0.000000) 

+13] = 1/2415602769592320(0.000000) 

= 32744383961/169885900800(0. 192743) [1520584461002805608448] 

-10] = 519961/2951427383938252800(0.000000) 

-9] = 254449/10540812085493760(0.000000) 

-8] = 544051261/491904563989708800(0.000000) 

-7] = 17186911411/885428215181475840(0.000000) 

-6] = 267372853813/1545985772539084800(0.000000) 

-5] = 9682910835401/3607300135924531200(0.000003) 

-4] = 10292911323827/885428215181475840(0.000012) 

-3] = 11916594797204821/97397103669962342400(0.000122) 

-2] = 84872899890985/61839430901563392(0.001372) 

-1] = 10174542900243355387/97397103669962342400(0.104465) 

+1] = 7257687613690421939/24349275917490585600(0.298066) 

+2] = 32905496396048207/32465701223320780800(0.001014) 

+3] = 294333194742383/4427141075907379200(0.000066) 

+4] = 630305608616617/32465701223320780800(0.000019) 

+5] = 378893651143/220855110362726400(0.000002) 
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+6] = 914315912029/2563081675525324800(0.000000) 
+7] = 1058275765039/32465701223320780800(0.000000) 
+8] = 10007111909/4058212652915097600(0.000000) 
+9] = 205655749/2563081675525324800(0.000000) 
+ 10] = 4511477/366154525075046400(0.000000) 
+11] = 2800937/5126163351050649600(0.000000) 
+ 12] = 33337/5126163351050649600(0.000000) 

= 943387771/18202060800(0.051829) [416456526393846005760] 

-10] = 33337/8321693751705600(0.000000) 

-9] = 83722897/295142738393825280(0.000000) 

-8] = 1149868206293/146095655504943513600(0.000000) 

-7] = 1532007576709/9130978469058969600(0.000000) 

-6] = 50123746848421/97397103669962342400(0.000001) 

-5] = 9301249107517/1739233994106470400(0.000005) 

-4] = 967852006786867/58438262201977405440(0.000017) 

-3] = 6951748079365363/73047827752471756800(0.000095) 

-2] = 36857804567275229/29219131100988702720(0.001261) 

-1] = 43645307985941263397/292191311009887027200(0.149372) 

+1] = 3660220040113447207/18261956938117939200(0.200429) 

+2] = 442368372384129727/292191311009887027200(0.001514) 

+3] = 489893102204383/2981543989896806400(0.000164) 

+4] = 7068099505148281/292191311009887027200(0.000024) 

+5] = 10499427847969/1014553163228774400(0.000010) 

+6] = 888556569941/869616997053235200(0.000001) 

+7] = 98029940871487/292191311009887027200(0.000000) 

+8] = 4599455666459/292191311009887027200(0.000000) 

+9] = 83722897/147571369196912640(0.000000) 

+10] = 33337/4160846875852800(0.000000) 

= 400879067/2316625920(0.173044) [1415737696661966684160] 
-6] = 6785814997/10821900407773593600(0.000000) 
-5] = 10357492793639/97397103669962342400(0.000000) 
-4] = 161442370258651/97397103669962342400(0.000002) 
-3] = 181501819795279/6493140244664156160(0.000028) 
-2] = 87969285921507127/97397103669962342400(0.000903) 
-1] = 66389099956276679/512616335105064960(0.129510) 
+1] = 2903737047336093503/9739710366996234240(0.298134) 
+2] = 13339984584863687/5410950203886796800(0.002465) 
+3] = 9234838162290949/24349275917490585600(0.000379) 
+4] = 88940575640699/1739233994106470400(0.000051) 
+5] = 10314251483003/1475713691969126400(0.000007) 
+6] = 75172645647961/48698551834981171200(0.000002) 
+7] = 9251203375/68348844680675328(0.000000) 
+8] = 15404457773/2705475101943398400(0.000000) 
+9] = 4004766797/32465701223320780800(0.000000) 
+10] = 903919/507276581614387200(0.000000) 
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V^56 
^56 

657 

V^i 

^57 
^57 
V^57 
V^57 

V^i 

V57 

^57 
^57 
^57 
V^57 
^57 
^57 
^57 
^57 

658 
^58 
^58 
V^58 
V^58 
^58 

V^s 

^58 
V^58 
^58 
^58 
^58 
^58 
^58 
V^58 
^58 
T^58 
^58 
V^58 
^58 
^58 



+11] = 403/26698767453388800(0.000000) 
+ 12] = 41/2029106326457548800(0.000000) 

= 30907469237/509657702400(0.060644) [505006549580537266176] 

-9] = 9899827187/292191311009887027200(0.000000) 
-8] = 2078227018577/292191311009887027200(0.000000) 
-7] = 279694033471/4869855183498117120(0.000000) 
-6] = 57975124388297/146095655504943513600(0.000000) 
-5] = 531949181968363/146095655504943513600(0.000004) 
-4] = 31272327932357/1987695993264537600(0.000016) 
-3] = 12585903870578537/97397103669962342400(0.000129) 
-2] = 57593752652897383/32465701223320780800(0.001774) 
-1] = 43587023408027517991/292191311009887027200(0.149173) 
+1] = 1530357192955689677/7304782775247175680(0.209501) 
+2] = 44651341400352107/24349275917490585600(0.001834) 
+3] = 777758638132831/4058212652915097600(0.000192) 
+4] = 148197040200331/10821900407773593600(0.000014) 
+5] = 1954397988395299/292191311009887027200(0.000007) 
+6] = 182673651588179/292191311009887027200(0.000001) 
+7] - 1351993586243/12174637958745292800(0.000000) 
+8] = 2078227018577/146095655504943513600(0.000000) 
+9] = 9899827187/146095655504943513600(0.000000) 

= 6267235837/70718054400(0.088623) [750950012179833421824] 
-9] = 324250979/2563081675525324800(0.000000) 

-8] = 23369547097/2563081675525324800(0.000000) 

-7] = 83894965921/512616335105064960(0.000000) 

-6] = 287011511657/732309050150092800(0.000000) 

-5] = 767530296899437/97397103669962342400(0.000008) 

-4] = 74756283949183/8854282151814758400(0.000008) 

-3] = 9067008775798753/97397103669962342400(0.000093) 

-2] = 122681701888511237/97397103669962342400(0.001260) 

-1] = 189447703516884708733/1412258003214453964800(0.134145) 

+1] = 207458630795906734361/941505335476302643200(0.220348) 

+2] = 18651227966840171/8116425305830195200(0.002298) 

+3] = 2862258848807693/16232850611660390400(0.000176) 

+4] = 573932994775373/16232850611660390400(0.000035) 

+5] = 11734282833971/2434927591749058560(0.000005) 

+6] = 26878790210113/48698551834981171200(0.000001) 

+7] = 23257660896641/97397103669962342400(0.000000) 

+8] = 215069764351/19479420733992468480(0.000000) 

+9] = 230825334239/48698551834981171200(0.000000) 

+10] = 5755065131/24349275917490585600(0.000000) 

+11] = 68623789/24349275917490585600(0.000000) 
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659 = 163742997419/3057946214400(0.053547) [27693220083749128765440] 



^59 


[-15] = 


= 41/6763687754858496000(0.000000) 


Vm 


[-14] = 


= 61/84546096935731200(0.000000) 




[-13] -- 


= 1319/48312055391846400(0.000000) 


1^59 


[-12] = 


= 2483197/6763687754858496000(0.000000) 


Vm 


[-11] = 


= 4760227907/1460956555049435136000(0.000000) 


Vm 


[-10] = 


= 36261195821/292191311009887027200(0.000000) 


Vm 


[-9] = 


99251424791/91309784690589696000(0.000000) 


1^59 


[-8] = 


145306469550229/5843826220197740544000(0.000000) 


V^59 


[-7] = 


96842410787587/922709403189116928000(0.000000) 


V^59 


[-6] = 


908924421985849/834832317171105792000(0.000001) 


Vm 


[-5] = 


15866293406797/2385235191917445120(0.000007) 


Vm 


[-4] = 


7538369215801399/273929354071769088000(0.000028) 


Vm 


[-3] = 


469300376053615651/1593770787326656512000(0.000294) 


Vm 


[-2] = 


30998833972906813663/17531478660593221632000(0.001768) 


Vm 


[-1] = 


101336333823239581902691/1034357240975000076288000(0.097970) 


Vm 


[+1] = 


15731441434639484641049/103435724097500007628800(0.152089) 


Vm 


[+2] = 


503906908616282993/318754157465331302400(0.001581) 


Vm 


[+3] = 


2029679255343050999/8765739330296610816000(0.000232) 


Vm 


[+4] = 


17982332246708579/626124237878329344000(0.000029) 


Vm 


[+5] = 


453773605688981/106251385821777100800(0.000004) 


Vm 


+6] = 


147689548546613/166966463434221158400(0.000001) 


Vm 


+7] = 


1221525888512777/17531478660593221632000(0.000000) 


Vm 


[+8] = 


116187185196623/5843826220197740544000(0.000000) 


Vm 


+9] = 


5560595677129/5843826220197740544000(0.000000) 


Vm 


+ 10] = 


= 147844102799/1460956555049435136000(0.000000) 


Vm 


+ 11] = 


= 1769351693/417416158585552896000(0.000000) 


Vm 


+ 12] = 


= 2483197/3381843877429248000(0.000000) 


Vm 


+ 13] = 


= 1319/24156027695923200(0.000000) 


Vm 


+14] . 


- 61/42273048467865600(0.000000) 


Vm 


+15] = 


= 41 /3381843877429248000(0.000000) 
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= 67410047209/254828851200(0.264531) [139128414846556255027200] 

-11] = 48572009/307569801063038976000(0.000000) 

-10] = 7564923041/102523267021012992000(0.000000) 

-9] = 241946098633/61513960212607795200(0.000000) 

-8] = 82710948197627/5843826220197740544000(0.000000) 

-7] = 618201928419193/1460956555049435136000(0.000000) 

-6] = 3617639494318727/1460956555049435136000(0.000002) 

-5] = 25317575928648679/1947942073399246848000(0.000013) 

-4] = 56398540189930703/531256929108885504000(0.000106) 

-3] = 93326831348910953/584382622019774054400(0.000160) 

-2] = 2254270194640917743/834832317171105792000(0.002700) 

-1] = 63503355003238436449/649314024466415616000(0.097801) 

+1] = 2119164534676506715771/5843826220197740544000(0.362633) 

+2] = 2258443180589756309/834832317171105792000(0.002705) 

+3] = 105935771427821093/531256929108885504000(0.000199) 

+4] = 540048225702181/39753919865290752000(0.000014) 

+5] = 3885941279688631/973971036699623424000(0.000004) 

+6] = 7952143409509/26088509911597056000(0.000000) 

+7] = 68408181635527/973971036699623424000(0.000000) 

+8] - 119925986148293/5843826220197740544000(0.000000) 

+9] = 5490125991767/5843826220197740544000(0.000000) 

+10] = 244847740997/307569801063038976000(0.000000) 

+11] = 14716614571/973971036699623424000(0.000000) 

+12] = 2237256377/243492759174905856000(0.000000) 

+13] = 30677699/177085643036295168000(0.000000) 
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= 19902302711/764486553600(0.026034) [863065094810581724037120] 

-13] = 33337/19223112566439936000(0.000000) 

-12] = 26683982221/18115861282612995686400(0.000000) 

-11] = 4775116095173/30193102137688326144000(0.000000) 

-10] = 69756035473169/15096551068844163072000(0.000000) 

-9] = 312279556434773/45289653206532489216000(0.000000) 

-8] = 72308884322139191/362317225652259913728000(0.000000) 

-7] = 34398025829829313/120772408550753304576000(0.000000) 

-6] = 457539002368786193/543475838478389870592000(0.000001) 

-5] = 2575014867714166963/181158612826129956864000(0.000014) 

-4] = 3066031235443749497/155278810993825677312000(0.000020) 

-3] = 15639813937655487379/67934479809798733824000(0.000230) 

-2] = 1423557203480883778099/1086951676956779741184000(0.001310) 

-1] = 2061939117042950979905/13958747851444960886784(0.147717) 

+1] = 11468365862330930545275961/66304052294363564212224000(0.172966) 

+2] = 44062736926197420397/24703447203563175936000(0.001784) 

+3] = 25978409273977745917/135868959619597467648000(0.000191) 

+4] = 19021961155258063699/1086951676956779741184000(0.000018) 

+5] = 89329015884779113/12077240855075330457600(0.000007) 

+6] - 2382219017953771/5200725727065931776000(0.000000) 

+7] = 4359398737953311/25879801832304279552000(0.000000) 

+8] = 185240586402109/1848557273736019968000(0.000000) 

+9] = 312280878311483/90579306413064978432000(0.000000) 

+10] = 69756035473169/30193102137688326144000(0.000000) 

+11] = 4775116095173/60386204275376652288000(0.000000) 

+12] = 26683982221/36231722565225991372800(0.000000) 

+13] = 33337/38446225132879872000(0.000000) 
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